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1.0  INTRODUCTION 

Space  based  military  optical  surveillance  systems  to  perform  a  variety  of 
missions  are  currently  under  consideration.  An  important  consideration  of  the 
design  of  such  a  system  is  the  choice  of  spectral  region.  Reliable  charac¬ 
terization  and  modeling  of  extremes  in  natural  backgrounds ,  and  of  khe  expected 
nuclear  perturbation,  is  necessary  to  optimally  choose  a  spectral  region. 

Here  we  report  an  effort  to  model  field  data  on  infrared  auroral  and  natural 
background  obtained  mainly  by  rocket  borne  sensors  in  the  A. 3,  2.7  and  LWIR 
spectral  regions,  to  predict  effects  on  structure,  and  to  conceive  field  experi¬ 
ments  to  most  effectively  enhance  the  data  base. 

Specifically  our  effort  was  concerned  with  development  of  a  two  dimensional 
time  dependent  radiation  transport  model  of  the  4.3pm  aurora.  Predictions 
of  spatial  structure  in  auroral  nadir  radiance,  and  experimental  concepts  for 
verification  were  addressed.  Nadir  radiance  gradients  in  the  CC^  4.3pm  aurora 
are  predicted  to  be  comparable  to  those  in  the  2.7pm  aurora  even  though  the 
DC  component  of  the  former  is  =10  times  the  brighter.  A  rocket  borne  arc 
over  flown  cryogenlcally  cooled  tilting  etalon  spectrometer  in  series  with  a 
CO^  absorption  cell  wheel  is  shown  to  be  the  optimum  experiment  to  provide 
verification.  The  modelling  facility  is  also  used  to  estimate  the  4.3pm  signal 
levels  that  might  be  observed  in  the  electron  gun  mother-  sensor  module 
daughter  spatial  EXCEDE  experiment  and  to  delineate  emission  mechanisms  for 
that  experiment.  A  detailed  investigation  of  the  CO,,  4.3pm  zenith  radiance 
data  obtained  under  undisturbed  night  time  conditions  on  4/11/74  revealed 
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compelling  evidence  for  the  existence  of  an  OH  airglow  related  mechanism 
OH(v)  +  N2  ♦  OH(v-l)  +  N*  followed  by  N*  +  C02  ♦  N2  +  C02(001)  followed  by 
partial  radiation  trapping  C02(001)  +■  CO^+  hv  4.3pm  for  producing  structured 
C02  4.3pm  radiance  under  undisturbed  night  time  conditions.  Detailed  window 
LWIR  C02  emission  earthlimb  synthetic  spectra  that  are  pumped  by  structural 
mechanisms  such  as  aurora,  airglow  and  fluorescence  are  also  presented.  The 
model  is  also  applied  to  very  high  sensitivity  auroral  4.3pm  zenith  radiance 
data  obtained  10/26/78.  These  data  are  shown  to  provide  evidence  for  a  prompt 
4.3pm  auroral  component  that  is  0.3%  energy  efficient  and  appears  to  be 
quenched  strongly  with  decreasing  altitude. 


2.0  THE  TWO  DIMENSIONAL  TIME  DEPENDENT  C02  4.3pm  AURORA,  RESULTS  AND  IMPLICATIONS 

2.1  Introduction 

Earlier  studies  (Kumer  1974,  1977a)  have  established  that  excitation  of  N2  by 
precipitating  particles  followed  by  N2  +  C02  ^  N2  +  C02(001)  followed  by  C02(001) 

-*C02  +  hv  (4.3pm)  is  the  principal  mechanism  for  generating  the  4.3pm  aurora. 

These  studies  have  also  shown  that  multiple  scattering  of  the  CO emission 
(C02 (001)-eC02  +  hv  (4.3pm))  significantly  impacts  the  temporal  and  altitude 
dependence  of  the  4.3pm  aurora.  Thus  we  world  also  expect  that  multiple 
scattering  would  tend  to  broaden  the  4.3pm  arc  with  respect  to  tht  3914A  arc. 

To  calculate  this  effect  it  is  necessary  to  solve  a  2  dimensional  time  depen¬ 
dent  radiation  transport  integral  equation.  We  have  done  this  in  application 
to  a  nominal  arc  model  in  which  the  arc  is  assumed  to  have  stabilized  in 
magnetic  latitude  for  2  minutes.  The  arc  is  assumed  to  have  been  moving  at 
3  km/min  in  a  direction  away  from  the  pole  prior  to  stabilization  time.  The 
e  fold  half  width  for  energy  deposition  is  assumed  to  be  ^*6.5  km.  Our  calcu¬ 
lation  shows  that  for  a  position  20  km  away  from  arc  center  in  the  direction 
of  motion  that  lateral  radiation  transport  increases  the  4.3pm  nadir  radiance 
by  approximately  a  factor  50  over  what  one  would  predict  if  lateral  transport 
were  neglected. 

Details  of  the  calculation  follow.  The  discussion  is  also  given  in  a  paper  by 
Kumer  (1978). 

2.2  Details  of  the  Calculation 

Two  spatial  coordinates  are  used  to  specify  the  arc  geometry,  namely,  the  alti¬ 
tude  z  and  the  distance  x  measured  in  the  north-south  direction  from  the 
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arc  center  (north  is  +) .  The  arc  is  assumed  to  extend  indefinitely  in  the 
y  direction  (East-West).  The  exact  2D  radiation  transport  equation  takes 
the  form 


S2D(z,  x)  -  SQ(z,  x) 


+  Jdx'dz'H^z,  z'  |x-x'|  )o[C02]z  n(z')S2D(zJ  x')  •  (2-1) 


The  quantities  S,  Sq,  Cl,  |C0oJ  and  E  are  defined  by  Kumer  and  James  (1974), 
which  will  be  referred  to  as  KJ.  The  2D  Green's  function  is  given  by 


where 


HoD  "  /V  H b(^)/%  **' 


•'  =  (x  -  x')‘  +  (z  -  i')c  +  y*‘ 


and  a  are  defined  in  KJ,  r'-  r  is  a  vector  from  the  point  r  =  (x,  0,  z) 
of  photon  absorption  to  the  point  r'  =  (x',  y' ,  z')  of  the  photon  emission, 
T)  is  the  column  density  of  CO^  molecules  along  the  vector  r  -  r'.  It  is 


useful  to  note  that 


=  H2D(°!  N  *  N'l’  I  X  ‘ 


where 


N  *  J,0°dz1[CO?]z  and  N'  -  dzJCQgl,  • 


Thus,  given  some  distribution  for  the  inhomogeneous  terms  Sc(z,x)  (the  initial 
excitation  rate),  we  solve  Eqn.  (2-1)  to  obtain  the  exact  two-dimensional 
solution  S2q(z,x) . 
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In  this  study  we  solved  the  problem  in  two  ways;  first  by  solving  the  exact, 
two-dimensiont 1  (2D)  integral  radiative  transport  equation  numerically. 

To  obtain  this  solution  by  iteration  requires  five  nested  loops  in  the 
computer.  We  also  solve  the  problem  by  introducing  an  approximate 
ID  radiative  transport  integral  equation  which  applies  for  each  altitude 


mesh  point  .  This  seccnd  metlod  requires  four  nested  loops  if  the 
ID  equations  are  solved  by  iteration  for  a  particlar  arc  geometry.  If  the 


inverse  G  =  (1  -  H  ) 
z  z 


is  obtained  at  each  z^ ,  then  3  nested  loops  are 


required  to  obtain  a  solution  for  any  particular  arc.  This  second  approach 


may  be  applied  tc  a  large  variety  of  arc  geometries. 


We  show  that  the  approximate  method  yields  results  that  are  very  close  to 
what  is  achieved  if  the  exact  method  is  employed.  For  many  purposes,  the 
accuracy  of  the  approximate  method  may  be  practically  sufficient  when  weighed 
against  (i)  the  extra  computer-time  cost  of  obtaining  an  exact  result  or 
(ii)  cases  where  inaccuracy  in  the  data  do  not  warrant  a  costly,  exact  solution. 

First  we  discuss  tie  iterative  technique  to  achieve  an  exact  solution  to 
Eqn.  (2-1).  The  computation  proceeds  most  rapidly  (and  economically)  if  we 
use  a  first  guess  trial  solution  (z, x)  that  is  as  realistic  as  possible. 

A 

An  excellent  trial  solution  S2D(z,x)  may  be  attained  in  a  2-step  process  as 
follows . 
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The  first  step  is  to  obtain  the  plane- parallel  solution  S^(z,x)  which  applies 
for  each  mesh  point  x  by  solving  the  1-dimensional  plane- parallel  transport 
equations 

SID(o  N,x)  -  S Q(o  N,*: 

+  J  CTdN'  Hhp(  |  N  -  N'|3)a(N')SU)(aN',x)  ,  (2-2) 

where  we  have  replaced  the  coordinate  a  by  N  =  j  dz*[CCh]  *  and  the  notation 

*J  £  Z 

z 

fo  l  Iowa  that  in  K.T.  The  do’ ails  of  lateral  transport  are  not 

accounted  for  in  executing  the  procedure  which  generates  S1D(oN,x);  therefore, 
the  x  dependence  will  be  incorrect.  However,  since  ^  *  j  dx  H^D» it  is  easy 

— <X> 

to  see  by  integrating  eq.  (l)  and  (2)  by  r°°  dx  that 

^-co 

J°°dx  Sld(oN,x)  =  J“  dx  S2D(on,x) 

-00  -00 

i j  an  identity;  due  to  lateral  transport, 

|  |jS1D(0li,x)|  >  I  ^S_,D(aN,x)| 

will  ho  id  rigorous Ly . 

The  next,  step  is  to  generate  an  appr oximotion  S^((TN,x)  that  is  the  solution 
of  an  approximate  l  dimensional  integral  transport  equation  for  each  mesh 
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•  -I*-*  • 


point  aN,  associated  with  z  .  That  is, 

J  J 

S2DaN,x)  '  So^N>x) 

+  n  (aN)  Jdx'H^D(aN,  |oc-x'  ')S^d(ctN,x')  ,  (2-3a) 


S^(aN,x)  =  (1  -  n(l-E))  S1D(aN,x)  , 

whore  hJd  is  given  by 

>4  ■  I  0"'  H2D 

and,  since  J  dx  H^D  =  1  -  E  the  identity 

J  <*  S2D  -  J*  dx  S1D  -  J  dx  s£d 

-ao  —00  —ao 

may  be  obtained  by  integrating  cq.  i.2-3)by  j  dx.  The  identity 

I  d»  S2D  -  I  «  4b 


(2-3b) 


(2-4) 


(2-3) 


i  :-h  i 


may  be  obtained  by  integrating  equs.  (2-1)  and  (2-3)  by  iadN  and  by  using  • .  n-.- 

J 

identity  JbdN  S,  =  J  adN  S ,^1-Q(  1-E) ) .  The  identities  (2-5)  and  (2-b)  art  useful 

A 

in  checking  on  the  accuracy  of  numerical  computations  of  S^,  S_,^  and  3.,^  .The 
analogous  equations  and  identities  fbr  the  3  "dimensional  case  are  given  in  Appendix  A. 

It  is  not  predictable  how  |-^  G^J  will  compare  with  |-j^  SgD|. 

In  Fig„(2-l)we  show  3  (z,x)/a  [C0„]^s  S^(oN,z)  which  is  utilized  describe  initial 
auroral  excitation  of  CO.  (OOl)  in  a  nominal  auroral  arc  as  described  in 
Appendix  B  below.  ..neso  :'uneticn.;  ’ire  given  for  too  •  highest  mesh  points. 
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z  —  92.5,  '/s .  1 ,  98.8,  107.2,  and  l4o  km.  The  full  sets  of  2 C  z  mesh  points 
and  31  x  mesh  points  are  given  in  Table  2-1. 
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Table  2-1 


<jN  AND  x  MET.1I  rOINTS  USED  IN  Tin-:  CALCULATION 
DESCRIBED  HEREIN 

(Note  that  .1400  +05  ■  0.14  x  10^) 


aN 


.1400+05 

.1897+00 

-100.0 

.  1072+05 

.8518+01 

-90.0 

.9884+02 

.4175+02 

-80.0 

.9508+02 

.9578+02 

-70.0 

.9252+02 

.1526+05 

-60.0 

.9096+02 

.2046+05 

-50.0 

.9018+02 

.2580+05 

-40.0 

.8318+02 

.5550+05 

-50.0 

.8584+02 

.6815+05 

-20.0 

.8022+02 

.1224+04 

-15.0 

.7168+02 

.1857+04 

-10.0 

.7605+02 

.2579+04 

-5.0 

.5105+0 


5000+02 
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are  shown  on  the  Interval  -30  >  ■  40  km 


LOG  to  (Sq) 


On  Fig.  2-2  we  plot  Sq(cjN,x)  vs.  x  for  z  =  95.1  and  107.2  km;  for  z  =  107.2  km , 
we  plot  s1D(aH,x),  S^d(oN,x)  and  S2D(aN,x). 

The  G2d  is  the  result  of  solving  eq.  2-1  by  iteration.  The  solution  0^ 
of  eq.  2-3  was  used  as  the  trial  solution.  Thirty  iterations  were  required 
to  achieve  convergence.  Convergence  was  accelerated  by  a  technique  described 
in  Appendix  C.  The  identity  2-5  is  satisfied  with  sufficient  accuracy 
(within  l/2#)  to  promote  confidence  that  there  is  no  major  logical  or  numerical 
error  in  the  computation.  The  integrations  in  eqs. (2-1),  (2-3)  and  (2-5)  were 
performed  over  the  range  -  100  <  x  <  240  km.  Outside  this  range,  we  assumed 
S(x^)  =  3(-100  km)  for  x  ^  -  100  km  and  S(x.  )  =  S(2U0  ion)  for  x^  S'  fto  km. 
Since  Cq(ctN,x)  is  a  maximum  at  x  *  0,  it  follows  that  the  integral  of  C 
over  the  range  x  =  -100  to  240  km  should  be  slightly  larger  than  the  integral 
of  SgD  or  on  the  same  interval.  This  we  found  to  be  the  case  by  about  1/27, 


From  Fig.  2-2  one 

sees  that, 

for 

z  = 

107.2 

km, 

ijL  0  I  .. 

>dx  °2D  "■ 

1- 

1  dx 

sn  I 

2D1 

(2-7) 

holds  for  all  x. 

The  "tail" 

in 

yiD 

for  x 

^  -  20  km  is  the  result 

of 

broadened  values 

which  occur 

in 

the 

nominal  arc  model  for  S^(aN,x) 

in 

the 

region  85  ^  z  =  97  km.  Essentially  the  same  kinds  of  results  are  plotted 


for  z  =  95.1  km  in  Fig.  2-3  Again  the  identity 2-5  is  satisfied  within 
0.67 For  z  =  95.1  km,  we  find  that 


N)u 


I  E  s2dI  <  I  e4>!- 


(2-8) 


In  fact,  the  inequality(2-8) is  true  for  all  z  mesh  points  except  z  =  107.2  km. 
The  single  qualitative  explanation  for  this  behaviour  is  (i)  that  is 
computed  on  the  assumption  that  ~  S(z,x)  =  0  and  (ii)  that  there  is  just 
one  maximum  in  the  quantities  So(z,x)  and  in  J  dx  Sq(z,x)  and  that  this 
maximum  occurs  at  z  =  107.2  km. 

2.3  Results 

On  Fig. 2-4  we  show  nadir  radiance  as  a  function  of  x  as  calculated  from  the 
various  solutions  S(aN,x)  cited  above  via 


M(x)  =  ar  dN'  T.  (aN')s(aN^x)  n  (lO  (2-9) 

D 

o 

00  A 

where  =  J  dz'fCOgl^/  and  =  50  km.  On  Fig.2-4it  is  seen  that  Rg^  =“  Rg^ 

(the  agreement  is  within  5$  or  better  on  the  region  of  x  shown  on  Fig. 2-4. 

Consideration  of  the  identity  2-6  shows  there  will  always  be  at  least  one  node 

A 

in  the  quantity  -  Sg^  on  the  interval  0  <  ah  <  oN^  and,  correspondingly, 

there  should  be  at  least  one  node  in  the  quantity  ul(S0^  -  Sg^)  on  this  interval. 

'This  would  tend  to  create  similarity  in  the  values  Rg^  and  .  for  ^he  4.5  Mm 

arc  model  that  we  consider  in  this  paper,  there  are  2  nodes  in  QT(£0^-S^) 

which  occur  on  either  side  of  the  maximum  in  Sg^.  This  should  further  increase 

the  similarity  of  RgD  and  Rg^.  Thus,  cue  might  indeed  expect  the  close 

A 

similarity  between  R^  and  Rg^  which  is  shown  on  Fig.  2-4. 


These  results  suggest  that  the  use  of  the  approximation  eq.  2-3  to  obtain 
approximate  plane  parallel  multidimensional  radiation  transport  solutions  may 
be  sufficiently  accurate  for  some  purposes. 
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NADIR  RADIANCE  IN  WATTS  /cm^sr 
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»rig.  2-4  The  Quantities  R^,  and  R^p  are  Shown. 
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2.4  Implications 

Lateral  gradients  in  the  nadir  radiance  R^Q  that  is  predicted  if  lateral 

transport  is  neglected  are  seen  to  be  considerably  steeper  for  all  values  of 

A 

x  than  for  the  predictions  or  Rjp  which  include  the  effects  of  lateral 
transport.  The  effect  of  lateral  transport  is  predicted  to  smoothout  lateral 
gradients.  This  effect  should  be  very  important  in  the  nuclear  case.  The 
effect  should  be  quantitatively  confirmed  experimentally  in  order  to  provide 
accurate  assessment  of  the  performance  of  advanced  nadir  viewing  scanner  sensors 
operating  against  the  auroral,  or  nuclear,  background. 


3.0  Experimental  Verification  of  the  Lateral  Photon  Transport  Smoothing 
Effect  in  the  CO^  4.3pm  Aurora 

3.1  Introduction 

In  this  chapter  we  provide  quantitative  analysis  for  several  alternate 
rocket  borne  experiments  to  verify  that  the  effect  of  lateral  radiation 
transport  of  C02  4.3  pm  photons  in  the  earths  optically  thick  upper  atmosphere 
renders  the  lateral  gradients  in  the  CO^  4.3  pm  auroral  nadir  radiance  no  more 
severe  than  the  lateral  gradients  in  auroral  2 . 7  pm  nadir  radiance  even  though 
the  DC  level  of  the  C02  4.3  pm  auroral  radiance  may  be  as  much  as  10  times 
that  of  the  2.7  pm  aurora  in  steady  state  conditions.  Verification  of  this 
point  is  important  for  assessing  high  altitude  nuclear  degradation  of  the  per¬ 
formance  of  advanced  nadir  viewing  scanner  systems  due  to  high  altitude  nuclear 
mechanisms  such  as  x-ray,  g-patch,  and  debris  patch  deposition  which  are 
similar  to  the  aurora. 

Persuing  this  goal,  we  considered  4  distinct  experimental  concepts  applicable 

for  the  4.3  pm  region,  and  briefly  considered  a  2.7  pm  experimental  concept  too. 

Most  of  these  concepts  involved  firing  a  rocket  over  an  auroral  arc  and  viewing 

the  aurora  in  the  nadir  in  order  to  measure  lateral  structure  in  auroral  nadir 

radiance  in  the  4.3  and  2.7  um  spectral  regions.  Nadir  viewing  4.3  pm  concepts 

included  a  Broad  Band  Filtered  Radiometer  (BBFR) ,  a  narrow  band  tilting  10 

transmission  channel  picket  fence  filtered  radiometer  (NBTFR)  and  the  use  of  HIRIS. 

In  each  case  the  instruments  are  to  be  used  with  a  wheel  of  C0„  absorption  cells 
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in  order  to  provide  the  added  high  resolution  ^  6  x  10  cm  near  the  centers  of 
the  C02  lines  that  is  necessary  to  discriminate  CO 2  4.3  pm  narrow  line  auroral 
nadir  radiance  from  the  earth  nadir  radiance  background.  It  turns  out  that  0.2  to 
0.3  cm  *  resolution  between  the  lines  is  also  necessary  in  order  to  suppress 
deleterious  effects  of  the  broad  band  LTE  thermal  stratospheric  contribution  to 
he  earth's  4,3  pm  nadir  radiance  background.  This  requirement  is  sufficient  to 
rule  out  the  BBFR  for  use  in  the  desired  arc  overfly  experiment.  The  NBTFR  and 
HIRIS  easily  satisfy  this  0.2  to  0.3  cm  1  spectral  requirement.  The  HIRIS  however 

turns  out  to  be  less  sensitive  than  is  required  for  a  successful  experiment. 

The  detailed  quantitative  analysis  presented  below  in  Sections  3.4  and  3.5  shows 
that  the  HIHTS  signal/noise  ^1  in  the  spatial  wing  region  of  the  nominal  arc 
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(approximately  15  km  or  more  south  of  arc  center).  However,  verification  of 
the  lateral  smoothing  effect  of  radiation  transport  on  the  CC>2  4.3pm  auroral 
arc  requires  accurate  measurement  in  the  spatial  wing  of  the  arc.  It  is  shown 
that  the  required  accuracy  may  indeed  be  provided  by  the  NBTFR  which  can 
achieve  signal/noise  of  approximately  30  in  the  spatial  wing  region  of  the 
nominal  arc. 

The  optimum  design  NBFTR  uses  liquid  Nitrogen  (U^)  cooled  optics  including 

the  tilting  filters,  and  LN-  cooled  baffling.  This  reduces  the  NEP  of  the 

.71  cm  diameter  InSb  detector  to  =3.0  x  10  w/t/Hz ,  nearly  the  inherent 

LN2  cooled  detector  noise  limit.  The  narrow  band  filter  is  a  4"  diameter 

Silicon  (index  of  refraction  =3.4)  substrated  Fabry  Perot  etalon  designed 

so  that  the  centers  of  10  of  its  transmission  channels  roughly  coincide  with 

12  16 

the  lines  p8,  plO,  ...,  p26  of  the  C  O2  001  -*•  000  band.  A  300  A  blocker 

filter  allows  just  the  10  channels  to  be  transmitted.  The  narrow  band  filter 

is  designed  so  that  the  transmission  channel  full  width  half  maximum 

FWHM  =  0.21  cm  ^  (i.e.,  4A)  for  an  FOV  with  the  1/2  cone  angle  a  =  2°.  The 

2 

through  put  of  the  optimum  NBTFR  is  =.062  cm  sr  for  an  overall  optics 

-13  2 

transmission  of  =  20%.  The  NER  =  5.0  x  10  w/cm  sr.  Spectral  scanning 

between  the  CO  line  centers  is  achieved  by  tilting  the  filter.  High  spectral 
1  -3-1 

resolution  of  =  6  x  10  cm  near  the  CO2  line  center  is  of  course  achieved 
by  utilization  of  the  CO2  absorption  cells.  Two  similar  NBTFR  instruments 
that  employ  cooled  optics  and  tilting  filter  have  been  developed  under 
DoD  contracts  DAAK40-75-C-0417  and  F3615-77-C-1190  for  deploument  in  the  Q 
hatch  of  the  U2. 

An  arc  underfly  verification  experiment  was  given  some  consideration.  Our 
impressions  on  this  concept  are:  1)  It  has  the  advantage  that  a  simple 
instrument  may  be  used  for  zenith  viewing  against  the  cold  skv  background 
however,  due  to  intervening  atmospheric  CO2  absorption  of  auroral  CO2  4.3um 
emission  it  is  probably  necessary  to  use  CO2  cells  in  order  to  accurately 
detect  weak  CO2  4.3pm  zenith  auroral  signals  from  rocket  altitudes  as  low  as 
85  km.  2)  It  has  three  disadvantages;  i)  the  rocket  spends  less  time  in 
the  altitude  range  where  useful  data  may  be  obtained  therefore,  decreasing  the 


23 


probability  that  phenomena  will  be  observed;  ii)  the  overfly  can  simul¬ 
taneously  monitor  energy  input  by  optical  measurements  and  by  measurements 
of  the  precipitating  auroral  particles;  iii)  the  overfly  measurements  are 
obtained  in  the  nadir  mode  which  is  exactly  similar  to  the  viewing  mode  of 
the  advanced  scanner  systems  that  the  data  will  support,  zenith  data  obtained 
in  the  underfly  experiment  will  not  share  this  advantage,  asthetic  as  it  may 
be.  As  part  of  our  study  we  also  verified  that  an  optimum  NBTFR  could  be  used 
effectively  in  an  arc  overfly  experiment  to  verify  that  lateral  gradients  in 
the  2.7pm  auroral  nadir  radiance  are  indeed  as  sharp  and  dramatic  as  what 
one  would  expect  from  an  essentially  prompt,  optically  thin  auroral  emission. 

3.2  Details  of  the  Requirement  to  Use  CO^  Absorption  Cells 

We  address  the  problem  of  devising  a  field  experiment  to  verify  that 
lateral  gradients  in  the  CO^  A. 3  pm  auroral  nadir  radiance  are  no  more  severe 
than  the  2.7  pm  nadir  radiance  even  though  the  DC  level  of  the  CC^  4.3  pm 
auroral  radiance  may  be  as  much  as  10  times  that  of  the  2.7  pm  aurora  in  steacy 
state  conditions.  Auroral  nadir  radiance  calculations  at  2.7  and  4.3  pm  for 
a  nominal  arc  as  described  in  Appendix  B  are  shown  in  Fig.  3-1.  The  lateral  gra¬ 
dients  in  the  auroral  CO^  4.3  pm  nadir  radiance  are  seen  to  actually  be  less  than 

in  the  auroral  2.7  pm  nadir  radiance  even  though  the  CO^  auroral  4.3  pm  nadir 
radiance  level  is  considerably  larger  than  at  2.7  ’.im.  The  nominal  30  kR  3914A 
aurora  moves  in  to  lateral  position  x  =  0  from  the  north  at  a  rate  3  km/min 
and  then  is  assumed  to  stabilize  for  2  minutes.  The  long  time  constant 
(  >  5  min)  for  auroral  CO  4.3  urn  emission  tends  to  smooth  the  4.3  Pm  lateral 

nadir  radiance  gradients  north  of  x  =  0.  This  long  time  dependence  for  the 

CC^  4.3  um  aurora  has  been  established  by  data  obtained  in  the  DNA/AFGL 
auroral  measurements  programs  and  by  the  detailed  analysis  of  these  data. 

The  smoothing  of  the  4.3  um  nadir  radiance  south  of  x  =  0  (i.e.,  for  x<0) 
is  predicted  to  be  the  result  of  lateral  radiation  transport.  This  effect 
has  not  as  yet  been  verified  experimentally.  The  work  described  in  this 
section  is  aimed  at  the  quantitative  analysis  of  various  experiments  which 
might  be  performed  in  order  to  verify  the  CO^  4.3  pm  lateral  nadir  radiance 
radiation  transport  smoothing  effect  which  is  depicted  on  Fig. 3-1.  Since  this 
effect  is  such  an  important  consideration  for  advanced  sensors,  it  is  clear 
that  some  kind  of  experimental  verification  is  necessary. 


The  verification  will  not  be  simple  in  an  arc  overfly  nadir  viewing  experiment 
due  to  the  earth  background.  This  point  is  illustrated  by  Fig.  3-2.  Auroral 
nadir  radiance  at  2.7  and  4.3  pm  appropriate  for  a  10  kR  3914A  steady  state 
auroral  arc  are  compared  with  a  night  time  1962  U.S.  Standard  Atmosphere  nadir 
radiance  background  calculation  (at  low  resolution  ~  0.01  urn)  on  Fig. 

The  background  levels  expected  to  obtain  at  2.7  and  4.3  pm  for  a  somewhat 

more  realistic  65°N  winter  atmosphere  (PKR  in  winter)  are  also  crudely 

drawn  in  as  dashed  lines  on  Fig  3-2.  Explicitly,  the  point  made  on  Fig.  3-2  is 

that  it  will  be  difficult  to  distinguish  the  auroral  nadir  radiance  from 

the  earth  background  nadir  radiance  if  resolution  like  0.01  Pm  is  used. 

It  is  also  clear  that  this  is  particularly  true  with  regard  to  nadir  measure¬ 
ment  of  the  4.3  pm  radiance. 

The  problem  has  a  potential  solution  if  the  nadir  radiance  is  viewed  in  high 
resolution.  This  point  is  illustrated  by  Fig.  3-3.  Fig.  3-3  shows  an  estimate 
for  the  detailed  spectral  structure  in  the  4.3  un  nadir  radiance  which  is 
calculated  for  auroral  conditions  similar  to  those  obtained  in  the  well  docu¬ 
mented  (AFGL-TR-77-0157 )  3/24/73  break  up. 

There  are  3  distinct  features  in  the  detailed  CC>2  4.3  pm  nadir  radiance 

spectrum,  auroral  R^,  stratopause  R^ ,  and  stratospheric  R^ . 

R.  :  ®  =  •  7 .8  :  40.  The  features  on  Fig.  3-3  are  easy  to 

A  p  S 

understand  qualitatively.  On  Fig.  3-4  we  show  the  CO^  v-j  vibrational  temperature 

altitude  profile  appropriate  for  the  3/24/73  event  and  on  Fig.  3-5  we  show 

the  altitude  z  where  optical  depth  =  1  for  a  given  displacement  in  cm  ^  from 

line  center.  It  is  apparent  that  R.  is  a  non-LTE  result  of  auroral  enhance- 

A 

ment  in  CO  2  \,^  vibration  temperature  which  occurs  in  the  altitude  region  z  > 

75  km  and  which  peaks  near  90  km,  that  R  is  due  to  the  LTE  v„  vibration 

P  3 

temperature  stratopause  maximum  which  occurs  near  the  z  =  50  km  region,  and 

that  Rg  results  from  pressure  broadened  emission  between  the  line  centers  in 

the  stratosphere  in  the  altitude  region  z  <  35  km.  The  detailed  spectrum  on 

-3  -1 

Fig.  3-3  shows  that  resolution  6  x  10  cm  or  better  is  required  to  dis¬ 
criminate  the  CO,  4.3  urn  aurora  from  the  stratopause  and  stratospheric  backgrounds. 
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Figure  3-4  The  CO.  v.  Temperature  Profile  Appropriate  For 
the  3/24/73  Aurora. 
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The  only  way  to  achieve  this  kind  of  resolution  with  state-of-the-art 
technology  involves  utilization  of  CO^  absorption  cells.  The  cells  may  be 
used  as  an  adjunct  resolution  element  by  positioning  them  in  front  of  a  broad 
band  (.2  pm)  filtered  radiometer,  (BBFR)  or  a  narrow  band  (4A  a  .21  cm-1)  picket 
fence  filter  tilting  radiometer,  (NBTFR)  or  a  HIRIS-llke  (resolution  *  0.5  cm”1 
9A)  interferometer. 

The  CO^  cells  may  be  used  to  discriminate  the  auroral  signal  as  follows.  He 

can  use  two  cells,  one  at  vacuum  and  another  with  an  amount  of  CO^  selected 

to  absorb  the  auroral  component  of  CO  4.3  pm  nadir  radiance  shown  on  Fig. 3-3. 

ij  1 

The  change  in  signal  R  on  interchange  of  the  two  cells  is  the  auroral 
nadir  radiance  component.  The  interchange  of  the  two  cells  allows  for 
weighting  the  auroral  C0^  4.3  pm  volume  emission  as  a  function  of  altitude 
as  is  currently  done  in  the  analogous  NASA  Nimbus  5  nadir  sounding  measurement 
of  15  pm  CO^  volume  emission  via  "The  Infrared  Temperature  Profile 
Radiometer  (1TPR)  Experiment"  which  is  described  by  Smith  et  al.  in  the  NASA 
Number  5  User's  Guide,  1972.  The  change  in  signal  R*^  on  cell  interchange  is 
given  by  an  integral  over  the  product  of  the  cell  interchange  weighting 
functions  W1^  and  the  CO^  volume  emission  Q  S  via 

R1"1  “  _1_  f  dz  wij  (z)  n  (*)  S  (z). 

4  V  J 

The  physics  details  of  the  calculation  of  the  cell  interchange  weighting  functions 
are  given  in  Appendix  G.  The  highlight  of  the  calculation  is  that  varia¬ 
tion  in  the  atmospheric  line  and  bandshapes  due  to  altitude  dependent  atmospheric 
temperature  variation,  and  pressure  variations,  are  accounted  for  in  detail. 


3.3  The  Broad  Band  Filtered  Radiometer  Concept 

Some  results  for  the  are  shown  on  Fig.  3-6  for  10  cm  length  cells  with 
various  CO^  and  broadener  gas  (Bf)  partial  pressures,  and  with  various  Lorentz 
line  broadening  parameters  as  listed  on  Fig.  3-6.  These  sounding  functions 
were  calculated  for  a  broad  band  (4.19  to  4.38  pm)  filtered  radiometer, 
but  they  can  also  be  used  to  gain  insight  as  to  requirements  of  resolution 
and  sensitivity  for  more  sophisticated  instrumentation  such  as  the  narrow 
band  picket  fence  tilting  filtered  radiometer,  or  a  HIRIS  like  instrument. 

12  1  12 

On  irig.  3-7  we  show  the  auroral  contribution  R  =7 —  [*dz  w  q  S  that 

A  *ir  *  a 

would  be  measured  in  an  auroral  arc  overfly  experiment  which  utilizes  a  4.19 
to  4.38  urn  broad  band  filtered  radiometer  (BBFR)  and  cells  1  and  2  from 
Fig.  3-6. 

On  first  consideration  it  might  seem  that  the  use  of  C02  absorption  cells  with 
the  BBFR  would  be  adequate  to  verify  the  lateral  transport  smoothing  effect  on 
CO 2  4.3  pm  auroral  nadir  radiance  in  an  arc  overfly  experiment.  There  are 
serious  problems  with  the  BBFR  arc  overfly  experiment  however.  These  problems 
are  caused  by  the  fact  that  the  broad  band  filter  does  not  spectrally  dis¬ 
criminate  against  the  large  stratospheric  nadir  radiance  background.  By 
inspection  of  Fig. 3-7  one  sees  that  for  x  s  -15  km  it  is  necessary  to  be  able 
to  accurately  discriminate  and  measure  the  auroral  contribution  to  nadir 

radiance  in  order  to  verify  the  CO,  4.3  um  lateral  transport  gradient  smooth- 

12 

ing  effect.  In  this  critical  region  the  auroral  contribution  R^  is 

12  -10  2 

quite  small  however,  R.  ^2.75  x  10  w/cm  sr  in  the  region  x  <  -  15  km. 

A  r* 

The  stratospheric  DC  contribution  R  to  the  nadir  radiance  observed  via  the 

s  -7  2 

BBFR  will  at  the  least  be  about  R  ~  2  x  10  w/cm  sr  (see  Fig.  3-2  above) , 

S  *" 

about  3  orders  of  magnitude  larger  than  the  auroral  nadir  radiance  we  wish  to 
detect . 

The  large  stratospheric  DC  signal  can  cause  problems  in  the  detection  of 
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Preliminary  Results  W  J  for  the  Sounding  Functions  Which  May  be  Achieved  by 
Interchange  of  Cells  Listed  Above. 
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Calculations  of  With,  and  Without, 

Including  the  Lateral  Radiation  Transport 
Smoothing  Effect. 


R,  in  two  ways*  First,  data  on  lateral  variations  in  the  stratospheric  com- 
ponent  of  4.3  pm  nadir  radiance  which  were  obtained  in  the  RM19  satellite  experi¬ 
ment  suggest  that  if  the  timing  of  cell  interchange  is  slow  (of  the  order  1/2  sec. 

or  slower)  the  lateral  structure  in  the  stratospheric  temperature  can  cause 

-2 

fluctuations  of  the  order  of  10  the  DC  signal  per  cell  interchange  and  this 

12 

will  be  more  than  enough  to  mask  the  auroral  signal  R  in  the  critical 
lateral  region  x  ^  -  15  km.  Secondly,  although  the  design  goal  would  be  to 
construct  the  cells  1  and  2  with  identical  transmission  properties,  it  may  well 
be  that  small  deviations  from  this  design  goal  could  cause  signal  changes  of 

-3 

the  order  10  of  the  stratospheric  background  on  cell  interchange.  Of  course 

preflight  calibration  of  the  transmission  properties  could  help  to  alleviate 

this  problem.  In  order  to  solve  both  of  these  problems  it  is  necessary  to  have, 

in  addition  to  the  C0„  absorption  cells  which  can  provide  a  resolution  capa- 

^  _  3  _i 

bility  of  a  few  times  10  cm  near  line  center,  an  additional  resolution  capability 
of  .2  to  .5  cm  1  (4  to  10  A)  between  the  lines  in  order  to  discriminate  against 
the  stratospheric  contribution  to  nadir  radiance  and  thereby  track  and  account  for 
the  lateral  fluctuations  in  the  stratospheric  component,  and  also  to  provide  in 
flight  calibration  of  the  transmission  characteristics  of  the  individual  cells. 

3.4  The  Narrow  Band  Solid  Etalon  Fabry  Perot  Tilting  Filtered  Radiometer  Concept 

With  this  purpose  in  mind  let  us  next  consider  the  use  of  CC>2  absorption  cells 
and  a  10  channel  narrow  band  solid  etalon  Fabry  Perot  tilting  filtered 
radiometer  (NBTFR) .  We  cite  a  similar  conceptual  instrument  described  by  Kumer 
and  Roche  (Arc  Overfly  Nadir  Radiance  Experiment  for  the  Multi  II,  Technical 
Brief,  LMSC  D457106,  May  1975).  A  sketch  of  the  design  for  an  improved 
version  of  the  NBTFR  for  use  in  the  arc  overfly  experiment  is  shown  on  Fig.  3-8  • 

This  instrument  is  considerably  improved  over  the  previous  conceptual 
design  which  is  given  in  the  1975  Kumer  and  Roche  technical  brief  in  that  the 
detector  and  all  of  the  optics  Including  the  narrow  band  tilting  filter  to  the 
right  hand  side  of  the  lens  LI  are  to  be  liquid  Nitrogen  (U^)  cooled. 

Similar  designs  for  LN^  cooled  optics  and  tilting  narrow  band  filters  have 
been  developed  under  DoD  support  in  two  background  measurements  programs, 
namely  HICAMP  (contract  No.  DAAK40-75-C-0417)  and  DBAR  (contract  No. 

F3615-77-C-1190) .  The  instrumentation  developed  under  these  programs  are 
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utilized  in  a  nadir  mode  from  high  flying  aircraft.  The  technology  for 
NBTFR  LN2  cooled  optics  including  tilting  filter  has  been  fabricated  and  is 
currently  in  the  debugging  process  under  the  auspices  of  these  programs. 

This  will  considerably  reduce  the  expense  and  the  risk  of  application  of 
similar  technology  to  the  arc  overfly  experiment. 

An  advantage  of  LN^  cooled  optics  including  tilting  narrow  band  filter  for  the 
arc  overfly  experiment  is  the  suppression  of  photon  noise  which  results  from 
nadir  viewing  of  the  earth  through  the  narrow  band  filter.  The  contribution 
to  NEP  due  to  the  250°K  earth  viewed  through  the  cold  narrow  band  filter, 
which  transmits  10  channels  of  width  0.2  cm  ^  each,  is  about  1.5  x  10  ^  w/*^Hz. 
The  inherent  InSb  0.71  cm  for  77°K  detector  noise  contribution  to  NEP  is  about 
2  x  10  w/v^Hz.  For  the  discussion  that  follows  we  use  NEP  =  3  x  10  ^ 
w/ /Hz  due  to  the  combined  effects  of  nadir  viewing  photon  noise  and  of  inherent 
detector  noise. 

In  the  Kumer  and  Roche  Technical  Brief  (1975)  it  was  shown  that  it  is  possible 

to  build  a  solid  Fabry  Perot  filter  with  a  10  channel  transmission  spectrum 

such  that  the  centers  of  the  transmission  channels  are  aligned  with  the  10 

12c  1602  001  -♦  000  band  lines  P8  through  P26  well  enough  so  that  the  net 

transmission  is  equivalent  to  7  completely  transmitted  lines.  To  achieve  this 

transmission  with  an  Si  substrate  (index  of  refraction  ~  3.4)  it  must  be  manu- 

3 

factured  with  a  thickness  tolerance  of  ±  1  part  in  10  ,  or  about  ±  1  urn.  In 
telephone  conversation  12-5-77  with  Perkin  Elmer  representative  Ed 
Strauss  Indicated  that  manufacture  of  an  Si  substrate  within  this 
tolerance  can  readily  be  achieved.  Use  of  the  Si  substrate  rather  than  quartz 
(index  of  refraction  ~  1.6)  allows  for  expanding  the  F0V  up  to  h  full  cone 
angle  a  -  2° .  With  Q  -  2°  the  full  width  half  maximum  (FWHMqo)  of  the  tr*”'-- 
mission  channels  is  broadened  by  a  factor  from  the  design  value  of  2  Jl  A 
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Rapid  filter  tilting  may  be  utilized  to  track  and  discriminate  against 

lateral  spatial  fluctuations  in  the  stratospheric  DC  signal  and  also  to 

provide  in  flight  calibration  of  the  transmission  characteristics  of  the  cells. 

Each  tilt  scan  should  be  progransned  so  that  about  2/3  of  the  tilt  scan  time 

period  t  is  used  to  displace  the  centers  of  the  transmission  spectra  from 

-4A  to  +4A  with  respect  to  the  band  line  centers,  thus  giving  an  effective 

line  center  match  up  dwell  time  of  about  1/3  t.  The  remaining  2/3  t  should 

be  used  to  obtain  a  rapid  tilt  scan  through  1/2  of  a  free  spectral  range  A-s 

in  order  to  detect  and  account  for  temporal  changes  in  the  stratospheric 

spectral  nadir  radiance  which  will  result  from  sensing  lateral  variation 

in  the  atmospheric  temperature  profile  as  the  result  of  the  lateral  rocket 

motion.  A  tilt  scan  period  of  about  1/10  second  can  readily  be  achieved. 

If  the  4  cells  listed  on  Figure  3-6  are  used,  with  interchange  each  1/2  sec. 

(take  0.1  sec  for  interchange  and  0.4  sec  dwell  time  for  each  cell)  then  the 

effective  on  line  dwell  time  for  any  pair  of  cells  during  a  complete  2 

second  cell  wheel  rotation  period  will  be  2  x  .4  x  (1/4)  x  (1/3)  sec  =  .067  sec. 

If  the  lateral  rocket  motion  is  ^  lkm/sec  then  the  corresponding  lateral 

spatial  resolution  is  ^  2km.  If  rocket  apogee  occurs  more  than  30km  above 

12 

the  centroid  altitude  (2  =“  105km)  for  W  one  might  consider  telescoping  the 
NBFTR  in  order  to  simultaneously  retain  both  the  o'  =  2°  l/2  cone  angle 
advantage  for  increased  throughout  and  the  2km  lateral  spatial  resolution 
(since  30  km  x  (4  ff/180)  2  km). 

Let  us  estimate  the  accuracy  to  which  we  can  obtain  the  CO2  auroral 
4.3  ^m  nadir  radiance  component  by  use  the  NBTFR  of  Fig.  3-8  sequenced  as  to 
tilt  scan  and  cell  interchange  as  described  above.  As  cited  above  the  10 
channel  picket  fence  filter  will  provide  100%  transmission  for  the  equivalent 
radiance  of  slightly  more  than  7  of  the  narrow  CO2  4.3  urn  auroral  nadir 

j 
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radiance  lines.  Information  provided  on  page  70  of  the  report  (AFCRL-TR- 

74-0334)  by  Kurner  indicates  that  there  are  about  40  strong  lines  in  the  auroral 

component  of  C02  4.3  p.m  nadir  radiance.  By  inspection  of  Fig.  3-7,  we  see 

12 

that  the  auroral  nadir  radiance  signal  in  a  single  line  on  interchange  of 

12  12 

cells  1  and  2  can  be  reasonably  estimated  by  /4G.  In  order  to  detect 

the  CO2  4.3  urn  nadir  radiance  in  the  crucial  lateral  transport  dominated 

12 

region  x  <  -15km  it  is  necessary  that  the  NBTFR  NER  <  <  7R^  in  that 

2  2 

region.  The  throughput  <£  of  the  NBTFR  is  given  by  *  Top  (tto?  )  (nr  ) 
where  Top  ^  20%  ,  a  is  the  FOV  1/2  cone  angle  2°,  and  r  is  the  aperture 

2  14  , — 

radius  2";  tp  =  .062  cm  sr.  The  NEP  is  approximately  3  x  10  w//Hz 
hence  the  NER  for  a  1  second  dwell  time  is  NER  —  NEP/y>^  5  x  10  ^ 
w/cm  sr/hz.  During  the  complete  2  second  cell  interchange  and  tilt 
scan  sequence  there  are  just  2  x  .067  seconds  of  on  line  dwell  time  available 
for  any  pair  of  cells  so  that  the  NER12  for  interchange  of  cells  1  and  2 

during  the  2  second  interchange  and  tile  scan  sequence  is  NER12  vj . 3  x  10~*2 

w/cm2  sr.  At  x  =  -15  km,  R*2  ^2.75  x  10_1°  w/cm2  sr  and  7R12/40  £4. 8  x 

1n-ll  .2  L 

iU  w/cm  sr-  Hence  in  the  crucial  region  x  <  -15  km  the  auroral  signal 

to  noise  Sa/N  which  can  be  obtained  by  use  of  the  NBTFR  sketched  on  Fig.  3-8 

along  with  the  4  C02  cells  listed  on  Fig.  3-6,  sequenced  as  to  tile  scan 

and  cell  interchange  in  2  second  intervals  as  described  above,  is  S  /N  < 

12  12  A  — 

7RL  /NER  ^37  for  x  £  -15  km.  By  inspection  of  Fig.  3-7  we  see  that 

Sa/N  retains  a  hefty  value  SA/N  242.6  at  x  =  -30  km.  Thus,  with 

2  km  lateral  resolution  in  the  region  -30  km  <  x  <  -15  km  we  can  obtain 

15  data  points  with  signal  to  noise  in  the  range  12.6  £  Sa/N  ^  37.  This  kind 

cf  Sa/N  in  the  critical  region  -30  £  x  £  -15  km  should  provide  for  unambiguous 

verification  of  the  smoothing  effect  of  lateral  radiation  transport  on  lateral 
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gradients  in  the  auroral  CO^  4.3ym  nadir  radiance.  It  also  becomes  clear 
that  a  footprint  as  large  as  8  km  will  not  degrade  the  verif ication,  hence 
„  the  NBFTR  need  not  necessarily  be  telescoped  in  order  to  obtain  useful  data, 
even  for  viewing  from  apogee  as  high  as  200  km. 


We  see  from  Fig.  3-7  that  nadir  radiance  due  to  the  mechanism  OH 
VV 
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N-  -  C09  (V,)  -*  C0_  +  h  "V.  will  begin  to  rival  that  of  auroral 

2  2  3  2  4.3  p.m 

12 

CO2  4.3  um  emission  in  contribution  to  the  measured  signal  for  x  in  the 
region  x  -30km.  The  use  of  cells  3  and  4  however, will  provide  enough 

12 

added  discrimination  against  the  OH  mechanism  so  that  its  contribution  R 

OH 

12 

to  the  measured  signal  R^  can  be  assessed  by  analysis  of  the  measured  signals 

14  34  12  12 

R^  and  R^  and  the  auroral  component  of  can  be  recovered  via  R^  - 

12  12 
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To  conlcude  our  discussion  of  the  use  of  the  NBTFR  along  with  CC^ 

absorption  cells  in  order  to  verify  the  lateral  radiation  transport  smoothing 

effect  on  lateral  gradients  in  the  CO2  4.3  pm  auroral  nadir  radiance  we  note 

12  12 

that  our  requirement  that  NER  <  <  7R^  is  satisfied  so  well  in  the 

critical  region  x  ^  -15km  that  we  can  in  fact  consider  relaxation  of  the 

cooled  black  baffling  requirements.  The  added  advantage  of  including  cells 

W  W 

3  and  4  for  discrimination  against  the  stratopause  and  the  OH(V)  N2  "* 

C09  (V„)  "•  CO-  4-  hV,  -  nadir  radiance  components,  as  well  as  the  filter 

tilting  discrimination  against  the  stratospheric  components  and  the  in  flight 
capability  for  filter  tilting  cell  transmission  characteristics  calibration, 
dictate  that  the  NBTFR  shown  on  Fig.  3-8  and  used  in  conjunction  with  the  CC^ 
absorption  cells  described  on  Fig.  3-6  is  eminently  feasible  for  verfication 
of  the  important  lateral  radiation  transport  emoothing  effect  on  lateral 
gradients  in  the  CC^  4.3  ^m  auroral  (and  therefore  for  the  similar  high 
altitude  nuclear  mechanisms)  nadir  radiance. 


41 


3.5  HIRIS  Concept 

Next  let  us  consider  the  use  of  HIRIS  in  conjunction  with  CO^  absorption  cells 
to  verify  the  lateral  radiation  transport  smoothing  effect  on  lateral 
gradients  in  auroral  CC^  4.3pm  nadir  radiance.  It  will  be  necessary  to 
include  a  4.21  to  4.32pm  LN^  cooled  filter  in  the  HIRIS  optics  in  order  to 
minimize  the  earth  shine  photon  fluctuation  contribution  to  HIRIS  NESR  and 
to  simultaneously  transmit  the  40  strong  auroral  nadir  radiance  lines. 
Resolution  of  0.5  cm  *  need  be  employed  to  discriminate  and  track  strato¬ 
spheric  nadir  radiance  fluctuations  and  to  provide  in-flight  calibration  of 
cell  transmission  characteristics.  For  the  purpose  of  estimating  S^/N  the 
auroral  signal  to  noise  that  may  be  achieved  in  the  critical  region 

x  <_  -15  km  by  using  the  HIRIS  we  use  the  HIRIS  parameters  listed  on  page  7 

of  the  report  AFGL-OP-TM-02  by  A.  T.  Stair,  J.  W.  Rogers  and  W.  R.  Williamson 
(1976),  and  we  get  NF.SR  ^  9  x  10  ^  w/cm^  sr  cm  1  for  HIRIS  viewing  the 

250°K  earth  through  an  LN  cooled  4.21  to  4.32pm  filter.  Details  of  the 

estimate  are  given  in  Appendix  H.  The  auroral  signal  per  CO^  auroral  nadir 
line  radiance  is  —  R^/.5  cm  *  ^  1 . 38  x  10  ^  w/cm^  sr  cm  *  at  x  =  -15  km. 
Thus,  even  if  all  40  lines  in  the  band  pass  are  considered,  the  HIRIS  signal 
to  noise  S.../N  at  the  critical  region  x  -15  km  is  Su./N  ^  1.  Hence, 
we  see  the  use  of  a  HIRIS  type  instrument  in  conjunction  with  CO^  cells 
would  not  provide  adequate  performance  to  verify  the  lateral  smooting 
effect  of  radiation  transport  in  the  spatial  wing  (x  -15  km)  of  the  CO2 
4.3pm  aurora. 
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The  telescoped  H1RIS  is  a  very  versatile  instrument.  For  example  it  could  be 
used  in  a  step  and  stare  nadir  viewing  mode  throughout  the  entire  Infrared 
spectral  region  in  a  long  term  satellite  borne  experiment  that  would  be 
extremely  useful  for  obtaining  comprehensive  high  spectral  resolution  AVD 
background  definition.  In  that  perspective  it  would  not  be  prudent  to  risk 
non-recovery  of  the  HIRIS  in  a  short  term  arc  overfly  experiment.  This  is 
especially  compelling  in  view  of  the  fact  that  the  highly  specific  NBTFR  is  a 
less  expensive  instrument (from  the  non-recovery  risk  point  of  view)  that  is 
considerably  more  sensitive  than  the  HIRIS  for  verification  of  the  lateral 
radiation  transport  smoothing  effects  on  the  lateral  CO^  4.3  ^m  auroral  nadir 
radiance  gradients . 

3.6  Arc  Underfly  Concept 

An  arc  underfly  experiment  might  be  useful  to  consider  as  an  alternate  approach 
for  providing  verification  of  the  lateral  radiation  transport  smoothing  effect 
on  lateral  gradients  in  CO^  4.3  ^m  auroral  nadir  radiance.  The  obvious 
advantage  is  that  viewing  the  aurora  in  the  zenith  sidesteps  the  nadir  viewing 
problem  of  discriminating  the  auroral  emission  from  the  earth  background.  How¬ 
ever,  now  the  aurora  must  be  viewed  through  large  optical  depths  unless  the 
sensor  is  at  90  to  95  km  altitude  or  higher.  On  the  other  hand,  once  the 
rocket  gets  too  high  (110  km)  the  auroral  signal  becomes  both  weak  and  diffi¬ 
cult  to  interpret.  Some  of  the  disadvantages  of  this  kind  of  experiment  are 
as  follows. 


Due  to  the  presence  of  the  0H(V)  mechanism  for  the  (X^  4.3  ^m  emission 

maximum  which  occurs  near  85  km,  and  to  the  large  optical  depth  of  atmospheric 
CO^  through  which  the  aurora  must  be  viewed,  it  appears  that  auroral 
zenith  signals  of  the  strength  found  in  the  criticai  region  x  <  -  15  km  will 
not  be  detectable  from  rocket  altitude  Z„  <  85  km,  even  if  CO,  cells  are  used 
to  attempt  to  sound  the  altitude  distribution  of  CO^  4.3  ^m  volume  emission 
above  the  rocket.  We  further  believe  that  data  obtained  from  rocket  altitude 
above  105  km  will  also  be  too  weak  and  uncertain  to  be  useful  in  the  x  s  -  15  km  j 

region.  Thus  useful  data  might  be  obtained  for  about  130  sec  as  the  rocket  i 

altitude  Zg  varies  from  85  to  apogee  (105  km)  and  then  back  to  85  km.  By 
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comparison,  in  an  overfly  NBTFR  experiment  the  rocket  can  be  programmed  to 
reach  apogee  at  220  km  and  will  obtain  useful  data  for  320  seconds.  Thus  in  the 
overfly  experiment  chances  for  observing  the  arc  in  the  critical  region  -40 
<  x  <  -15  km  are  considerably  enhanced  over  those  of  the  arc  underfly  experiment. 
Another  advantage  of  the  arc  overfly  experiment  is  that  energy  deposition  can 
be  monitored  both  optically  (by  nadir  measurements  in  the  Nj  LBH  system  for 
example)  and  also  by  observing  the  precipitating  auroral  particles.  Finally, 
the  nadir  viewing  arc  overfly  concept  is  more  appealing  esthetically  since  the 
advanced  earth  disc  scanner  systems  for  which  the  data  will  provide  support, 
will  also  be  used  in  a  nadir  viewing  mode. 

Thus  we  conclude  that  the  arc  underfly  exepriment  appears  to  be  considerably 
less  desirable  than  an  arc  overfly  experiment  which  utilizes  the  NBTFR  as 
described  above  for  the  purpose  of  verification  of  the  smoothing  effect  in  the 
lateral  CO^  4.3  pm  nadir  radiance  gradient  which  results  from  lateral  transport. 

3.7  An  NBTFR  2.7  pm  Arc  Overfly  Concept 

Finally,  we  should  address  the  verification  of  relatively  sharp  lateral 
gradients  in  auroral  2.7  p-m  nadir  radiance.  The  2.7  instrument  discussed 
by  Kumer  and  Roche  (1975)  can  be  improved  in  NER  by  a  factor  x  1/6  by  use  of 

the  improved  NBTFR  technology  discussed  above, namely  cooled  optics  and 

tilting  filter  and  the  use  of  an  Si  substrate  for  the  narrow  band  filter.  On 

the  other  hand  the  2.7  pra  auroral  signal  in  the  Kumer  and  Roche  Technical 

Brief  (1975)  is  over  estimated  by  a  factor  3  since  an  energy  efficiency  27. 
was  used  rather  than  0.7%  which  now  appears  more  realistic  in  view  of  recent 
independent  determinations  frcxn  a  large  number  of  sources  as  reported  at  the 
April  1979  DNA/AFGL,  HAES  Infrared  Data  Review  Meeting.  The  net  effect 
is  that  sensitivity  to  the  2.7  pa  auroral  nadir  radiance  that  is  a  factor 
x  2  better  than  depicted  in  the  Kumer  and  Roche  article  can  be  achieved  with 
a  NBTFR.  Performance  with  this  sensitivity  for  an  arc  similar  to  that  over¬ 
flown  on  3/27/73  in  the  Ice  Cap  program  is  illustrated  on  Fig.  3-9. 
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Figure  3- 


9  Illustrates  Anticipated  Performance  of  an  Optimized  2.8pm  NBTFR 
on  Overflight  of  an  Arc  Modelled  After  the  3/27/73  Arc 
Encountered  in  the  ICECAP  Program. 


3.8  Conclusions 


We  have  considered  4  distinct  experimental  concepts  to  verify  that  the  effect 
of  lateral  radiation  transport  of  COj  4.3  pm  photons  in  the  earths  optically  thick 
upper  atmosphere  renders  the  lateral  gradients  in  the  CO^  4.3  pm  auroral  nadir 
radiance  no  more  severe  than  the  lateral  gradients  in  auroral  2.7  pm  nadir 
radiance  even  though  the  DC  level  of  the  CC^  4.3  pm  auroral  radiance  may  be  as 
much  as  10  times  that  of  the  2.7  pm  aurora  in  steady  state  conditions.  Most  of 
these  concepts  involved  firing  a  rocketborne  sensor  over  an  auroral  arc  and  viewing 
the  aurora  in  the  nadir  in  order  to  measure  lateral  structure  in  auroral  nadir 


radiance  in  the  4.3  and  2.7  .m  spectral  regions.  Nadir  viewing  4.3  pm  concepts 
included  a  Broad  Band  Filtered  Radiometer  (BBFR) ,  a  narrow  band  tilting  10 
transmission  channel  picket  fence  filtered  radiometer  (NBTFR)  and  the  use  of 
HIRIS.  In  each  case  the  instruments  are  to  be  used  with  a  wheel  of  CO 

-3  -l 

absorption  cells  in  order  to  provide  the  added  high  resolution  <  6  x  10  cm 
near  the  centers  of  the  CO  2  lines  that  is  necessary  to  discriminate  COj  4.3  pm 
narrow  line  auroral  nadir  radiance  from  the  earth  nadir  radiance  background. 

It  is  also  necessary  to  have  0.2  to  0.5  cm  *  resolution  between  the  lines  in 
order  to  suppress  deleterious  effects  of  the  broad  band  LTE  thermal  stratospheric 
contribution  to  the  earth's  4.3  pm  nadir  radiance  background.  This  requirement 
is  sufficient  to  rule  out  the  BBFR  for  use  in  the  desired  arc  overfly  experiment. 
The  NBTFR  and  HIRIS  easily  satisfy  this  0.2  to  0.5  cm  1  spectral  requirement. 


The  HIRIS  however  turns  out  to  be  significantly  less  sensitive  and  more  costly 
than  the  optimum  design  NBTFR. 


The  optimum  design  NBFTR  uses  liquid  Nitrogen  (L^)  cooled  optics  including 

the  tilting  filters,  and  LN_  cooled  baffling.  This  reduces  the  NEP  of  the 

^  *“  1 A 

.71  cm  diameter  in  Sb  detector  to  =  3  x  10  w/ /  Hz,  nearly  the  inherent 
LNj  cooled  TnSb  detector  noise  limit.  The  narrow  band  filter  is  a  4"  diameter 
Silicon  (index  of  refraction  =  3.4)  substrated  Fabry  Perot  etalon  designed  so 
that  the  centers  of  10  of  its  transmission  channels  roughly  coincide  with  the 
lines  p8,  plO,  ...,  p26  of  the  ^C^O^  001  -*•  000  band.  A  300  A  blocker  filter 
allows  just  the  10  channels  to  be  transmitted.  The  narrow  band  filter  is 


46 


designed  so  that  the  transmission  channel  full  width  half  maximum  FWHM  0.21  cm 

(i.e.,  4A)  for  an  FOV  with  the  1/2  cone  angle  -  2°.  The  through  put  of  the 

2 

optimum  NBTFR  is  .031  cm  sr  for  an  overall  optica  transmission  of  102. 

-13  2 

The  NER  —  7.10  x  10  w/cra  sr.  Spectral  scanning  between  the  CO  line  centers 

1  -3  -1 

is  achieved  by  tilting  the  filter.  High  spectral  resolution  of  6  x  10  cm 
near  the  CO^  line  center  is  achieved  by  utilization  of  the  CO^  absorption  cells. 

k 

An  arc  underfly  verification  experiment  was  also  considered  and  rejected 
on  the  basis  that;  i)  the  rocket  spends  less  ti*me  in  the  altitude  range  where 
useful  data  may  be  obtained  therefore  decreasing  the  probability  that  phenomena 
will  be  observed;  ii)  the  overfly  can  simultaneously  monitor  energy  input  by 
optical  measurements  and  by  measurements  of  the  precipitating  auroral  particles; 
iii)  the  overfly  measurements  are  obtained  in  the  nadir  mode  which  is  exactly 
similar  to  the  viewing  mode  of  the  advanced  scanner  sensors  that  the  data  will 
support,  zenith  data  obtained  in  the  underfly  experiment  will  not  share  this 
advantage. 

As  part  of  our  study  we  also  verified  that  an  optimum  NBTFR  could  be  used 
effectively  in  an  arc  overfly  experiment  to  verify  that  lateral  gradients  in 
the  2.7  m  auroral  nadir  radiance  are  indeed  as  sharp  and  dramatic  as 
one  would  expect  from  an  essentially  prompt , optically  thin  auroral  emission. 
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4.0  APPLICATION  TO  PLANNING  FOR  SPATIAL  EXCEDE 

In  this  chapter  we  generate  estimates  for  the  4.3  ym  signal  levels  that  we 
expect  to  see  in  the  EXCEDE  mother-daughter  spatial  experiment  as  described  to  us 
by  Bob  O'Neil.  In  this  experiment  the  mother  vehicle  will  carry  the  electron  gun 
and  provide  at  altitude  z  -  90  km  a  120  kw  dosed  region.  Observations  will  be 
performed  from  the  nearby  daughter  vehicle.  For  the  purpose  of  calculation,  the 
dosed  region  can  be  approximated  by  a  magnetic  field  aligned  rectangle  of  70  m 
length  L  parallel  to  the  magnetic  field  and  a  20  m  width  D  which  is  perpendicular 
to  the  magnetic  field  but  parallel  to  the  horizontal  component  of  rocket  velocity. 

The  rectangle  is  assumed  to  have  a  one  cm  depth.  The  rectangle  has  velocity  com¬ 
ponent  V  =  600  m/sec  in  the  direction  perpendicular  to  the  magnetic  field. 

On  fig. 4-1  we  show  a  schematic  of  the  instantaneous  deposition  region  (IDR)  and 
adjacent  regions  as  viewed  from  the  daughter  vehicle  which  is  presumed  to  be 
moving  parallel  to  the  EXCEDE  electron  gun.  Adjacent  regions  are  the  £redosed 
region  (PR),  the  advanced  transport  region  (ATR)  and  the  _lower  transport  region 
(LTR).  We  shall  see  below  that  transport  rather  than  predosing  dominates  4.3  um 
radiance  due  to  prompt  electron  beam  induced  direct  excitation  of  CO^  in  the  state. 
We  call  this  4.3  pm  excitation  mechanism  I.  In  discussion  of  mechanism  I  we  refer 
to  the  predosed  region  (PR)  as  the  predosed  transport  region  (PDR) . 

We  expect  4  obvious  components  of  4.5  pm  radiation  to  be  present.  The  first  mech¬ 
anism  (Mechanism  I)  is  prompt  electron  beam  induced  direct  excitation  of  CO^  (v^) . 
Mechanism  I  will  produce  a  component  that  will  radiatively  relax  at  A»  419/sec, 
the  deposition  rectangle  will  move  one  D  perpendicular  to  the  field  line  with  a 
frequency  v/D  «  30/sec  so  this  component  will  be  essentially  prompt.  A  2  cd 
component  R**  will  result  from  prompt  electron  beam  induced  excitation  of  N* 
followed  by  the  reactions 
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Figure  4-1  A  Schematic  of  the  Instantaneous  Deposition  Region  and 
Adjacent  Regions  as  Seen  From  the  Daughter  Vehicle. 
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(4-n 


N*  +  C0„  -»  N.  +  CO^(OOl) 

c  tC  L- 

CO  ,(001)  -*  C0„  +  hv.  (4-2) 

This  second  mechanism  (Mechanism  II)  wilL  leave  a  long  DC  tail  behind  the 
deposition  region  because  the  time  constant  for  reaction  (4-1)  is  of  the  order 
of  a  minute  at  90  km  altitude. 


Radiation  transport 

hv  +  C0o  -*  C0n(001) 

will  also  cause  a  glow  outside  of  the  dosed  and  predosed  regions.  Both 
mechanisms  I  and  II  will  contribute  to  the  glow  in  the  transport  region  but, 
as  we  shall  see  below,  mechanism  II  will  dominate. 


We  also  considered  another  component  r^q+,  namely  4.}  urn  emission  from  N0+ 
produced  in  vibrational ly  excited  states.  At  90  km  altitude  this  N0+  emission 
will  be  heavily  quenched  by  so  that  it  will  be  essentially  prompt.  To 
generate  the  contribution  Rjj0+  from  N0+(v)  -*  N0+(v-i)  +  hv  we  followed  Tarr 
and  Archer  (DNA  }C97F)  in  assuming  a  production  eNQ+.  of  5  photons  per  N0+ 
produced,  we  used  rates  kj(v)  for 


Np  +  N0+(v)k}iV'1i^  +  N0+(v-l) 


(4-3) 


inferred  from  the  data  for  k^(l)  given  by  Gordietz  et.  al.  (1978)  and 
from  calculations  of  kj(v)/k^(l)  as  listed  in 

DNA5297F,  and  we  used  Einstein  coefficients  A^  ^  for  N0+(v)  -*  NQ+(v-l)  +  hv 
taken  from  Frank  Billingsleys  publications.  The  results  we  generated  on  this 
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bus  Is  Indicated  that  R^+  would  dominate  heavily  In  the  region  IDR.  However, 
Herb  Mitchell  reports  (private  communication)  that  in  earlier  electron  gun 
experiments  the  ratio  of  1*.}  pm  to  591  remains  constant  over  the  altitude 
region  80  •  z  ••  LO  km.  N*  W4A  emission  should  not  be  quenched  in  this 
region  but  R^+  should  be  heavily  quenched,  i.e.,  A^  ^/(A^  ^  +  ^(v)[N?l)  S  l/c 
for  z  '•  l.'O  km.  Thus  in  order  for  the  constant  ratio  591!+aA*3  dm  to  obtain 
in  the  region  30  <  z  <  l.?0  km  it  must  be  that  R^0+  is 

not  dominant..  For  this  to  occur  would  mean  that  either  ®Ng+  <  k^(l)  Is 

greater  than  given  by  Gordie tz  et.  a].,  A^  is  less  than  calculated  by 
Billingsley,  or  any  combination  of  the  above. 


For  our  purpose  ve  have  generated  estimates  RN0+  the  data  base  cited  above 
(eNg+  and  k.  (v)  k.(l)  from  PNA  5:H>7P,  k  ^  ( 1 )  from  Gordietz  et  .  al. ,  and  A^  ^ 
from  Hi  1 i ings ley) .  The  EXCEDE  mother -daughter  experiment  may  prove  these 
estimates  to  be  cons Iderab 1 y  erroneous.  This  will  indicate  deficiencies  in 
the  data  base,  rather  than  in  the  met  hods  we  use  to  estimate  R„,.  from  the 
data  base. 


Finally  there  is  the  ambient,  background  CO,  ^.5  uni  emission  to  consider.  At 
"Xt  km  this  will  be  dominated  by  the  Otf*  V  VcO„(OOl)  -»  CO,  +  hv  mechanism. 

Then  it  can  be  quite  variable.  This  mechanism  for  high  altitude  CO,  emission 

is  discussed  in  a  paper  by  Kumer  et  al.  ( 197S) . 


;e  suits  are  summarized  in  Table  1.  Clearly  R^  +  dominates  the  region  IDR 


and  R^  dominates  the  transport  regions  and  the  region  PR.  From  the  data  in 

IDR  we  can  get  ,  since  e  .  “  R.t/_+,  from  data  in  the  region  PR  we  can  get 
N0+  NO  NO  0 

C ^  ,  the  efficiency  for  prompt  high  energetic  election  collision  induced 

excitation  of  N*  by  means  other  than  vv  transfer  from  NO+(v)  since 
II 

€  +  €ma+  “  R  .  Finally  we  may  get  some  information  on  e  ,  the  prompt 

Ww  COp 

high  energetic  electron  collision  induced  efficiency  for  producing  COo(v^)  via 
e  +  C0.N  -»  e  +  CO,(v^)  from  data  in  the  transport  region  since  in  these  regions 
R^/R*^  ~  0. ;':V  based  on  the  efficiencies  we  have  assumed  in  order  to  generate 


these  preliminary  estimates. 


Details  of  the  calculation  of  the  estimates  in  table  4.1  are  given  in  Appen¬ 
dix  D. 
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5.0  EVIDENCE  FOR  AN  OH  V  N2  -*■  CO  (v3>-*-  C02  +  hv  (4.3pm)  Mechanism  for 
4.3pm  Airglow 

5.1  Introduction 

In  the  previous  work  by  Kumer  and  James  (1974)  and  Kumer  (1974)  it  had  been 
assumed  that  the  only  important  undisturbed  night  time  mechanisms  in  the 
altitude  region  z  >  50  km  for  producing  C02  excited  in  the  normal  mode  of 
vibration  were  (i)  excitation  by  thermal  collisions  and  (ii)  excitation  by 
the  radiative  absorption  of  4.3pm  earthshine  by  C02-  The  radiative  decay  of 
CO^tVj)  results  in  emission  near  4.3pm.  Therefore,  measurements  of  the 
undisturbed  night  time  atmospheric  4.3pm  zenith  radiance  for  altitudes  z  >50  km 
might  be  used  to  test  if  these  are  indeed  the  only  important  mechanisms  for 
producing  CO^fv^)  in  undisturbed  night  time  conditions. 

Since  1974  preliminary  4.3pm  zenith  radiance  data  sets  obtained  27  April  1973, 
25  February  1974,  11  April  1974,  6  March  1975,  and  12  March  1975  have  shown 
evidence  (Kumer,  HAES  Report  No.  57,  1976)  for  a  non-auroral  feature  near  85  km 
which  cannot  easily  be  explained  on  the  basis  of  mechanisms  (i)  and  (ii).  The 
consistent  appearance  of  this  feature  in  the  4.3pm  zenith  radiance  data  near 
85  km,  and  the  large  values  for  hte  rate  constants  (v)  (i  =  1  for  N2,  i  =  2 
for  02)  for  the  reactions 

0«(v)  +g)i“V>  OH(v-l)  + 

which  have  been  reported  by  Streit  and  Johnston  (1976),  had  stimulated  a 
preliminary  study  (Kumer,  1976)  to  investigate  the  mechanism 


OH(v)  +  N2 


OH(v-l) 


(5-2) 


followed  by  near  resonant  w  transfer  back  and  forth  between  N2  and  C02 


n  +  co2  t  n2  +  co2(ooi) 


and  4.3  pm  radiation  trapping 

C02(001)q*  C02  +  h“V 

as  the  source  of  the  feature  which  is  observed  near  80  to  85  km  altitude.  In  this 
chapter  we  show  in  detail  how  this  mechanism  explains  the  presence  of  the  feature 
which  occurs  near  85  km  in  the  data  set  obtained  on  11  April  1974.  The  work  this 
chapter  covers  has  been  published  in  J.  Geophys.  Res.  (Kumer,  et  al,  1978). 

5.2  ANALYSIS 

One  of  the  goals  for  the  1974  auroral  measurements  program  was  to  obtain 
a  4.3  pm  zenith  radiance  profile  under  quiet  conditions  so  that  there  would 
be  no  significant  auroral  component  in  the  data  thus  obtained.  The  evening 
of  11  April  1974  was  free  of  auroral  activity  so  the  "quiet  conditions" 
rocket  was  launched  that  evening.  The  4.3  pm  zenith  radiance  data 
are  shown  by  curve  D  on  Figure  5-1.  These  data  were  obtained  with  a  rocket- 
borne  liquid  nitrogen  cooled  CVF  spectrometer,  the  description  of  which  has 
been  given  by  Stair  et  al.,  (1973)  and  Wyatt  and  Frodsham  (1977).  Spectra  thus 
obtained  showed  an  emission  feature  near  4.3  pm  which  may  be  identified  as  C02Vj 
band  emission.  A  preliminary  version  of  these  data  were  reported  by  Stair  et  al. 
(1975)  and  the  final  version  of  these  data  are  reported  by  Wheeler  et  al.,  (1976). 

Data  above  92  km  are  lost  since  the  signal  is  smaller  than  instrument  noise 
for  z  >  92  km. 

The  curve  A  which  is  also  shown  on  Figure 5-lis  a  calculation  which 
accounts  for  just  the  two  traditional  nighttime  mechanisms  (Kumer  and  James, 

1974)  for  the  excitation  of  C02  (v.j) ,  namely  (i)  vibrational  excitation  by 
thermal  collisions  and  (ii)  absorption  of  4.3  pm  earthshine  by  C02<  In  order 
to  perform  this  calculation  it  is  necessary  that  we  use  the  most  accurate 
model  for  the  atmospheric  temperature  and  pressure  on  the  evening  of  ll  April 
1974  that  can  be  constructed  on  the  basis  of  whatever  pertinent  data  are 


v*  ■ 
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Fig.  5-1.  The  4/11/74  4.3  urn  zenith  radiance  data  are  shown  by  curve  D. 

Error  bars  Indicate  the  uncertainty.  Curve  A  shows  a  4.3  pm  zenith  radiance 
calculation  that  is  based  on  the  traditional  night-time  C02(v3)  excitation 
mechanisms  (i)  and  (ii)  which  are  discussed  in  the  text,  and  on  our  nominal 
atmospheric  model  which  is  also  described  in  the  text. 
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available.  Up  Co  60  km  altitude  our  model  is  based  on  Che  U.S.  Army  ACmospheric 
Sciences  Laboracory  meteorological  rockeC  data  which  were  obtained  in  a  sounding 
rocket  experiment  launched  from  the  Poker  Flat  Range  (PKR)  on  the  evening  of 
11  April  1974.  These  data  were  provided  by  R.  0.  Olsen  (1976).  Above  60  km  we 
joined  this  data  smoothly  to  a  data  point  at  87  km  which  is  interpolated  from 
OH  airglow  rotational  temperature  determined  (Ware,  private  communication  1978) 
from  interferometric  data  obtained  on  4/9/84  and  4/13/74  via  ground  based  instru¬ 
mentation  located  near  PKR.  No  such  data  were  available  from  the  10th  to  the  12th 
due  to  clouds  and  snowstorm  activity.  Thus  we  interpolate  the  mesopause  temperature 
T  -  203  +25°K  to  apply  for  4/11/74  on  the  basis  of  these  preliminary  data  supplied 

by  Ware: 


time(UT) 

date 

OH  rotational 
temperature 

1:30 

4/9/74 

205.9 

8:33 

4/9/74 

194.0 

11:00 

4/13/74 

207.6 

The  large  uncertainty  we  assign  to  the  interpolated  mesopause  temperature 
data  point  accounts  for  large  daily  and  hourly  variations  (Ware  et  al.,  1976;  and 
Noxon,  1978)  in  airglow  rotational  temperature,  and  presumably  in  mesopause 
temperature  as  well.  The  temperature  profile  TN(z)  for  our  nominal  model 
atmosphere  which  is  constructed  on  the  basis  of  available  data  as  explained 
above  is  shown  on  Figure  5-2. 

We  use  the  nominal  temperature  model  T^(z)  and  the  mechanisms  (i)  thermal 
collisions  and  (ii)  4.3  pm  earthshine  absorption  in  order  to  generate  the 
calculation  of  4.3  m  zenith  radiance  which  is  shown  by  curve  A  on  Figure  5-1. 
The  calculation  and  data  (curve  D  on  Figure  5-1  are  within  a  factor  0.78  agree¬ 
ment  below  z  v  60  km  altitude.  Considering  that  the  overall  calibration  of  the 
CVF  is  about  ±  25%,  and  also  considering  that  an  overall  change  of  2.5°K  in 
the  temperature  profile  can  cause  a  10%  change  in  4.3  pm  zenith  radiance,  the 
agreement  for  z  ^  60  km  between  the  data  and  calculations  as  shown  on  Figure  5-1 
is  quite  good.  Centered  at  about  85  km  there  seems  to  be  a  contribution  of 
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approximately  4  x  10  w  cm  *sr  4.3  .wn  zenith  radiance  data  that  can  not  be 
explained  by  our  nominal  model  atmosphere  and  the  traditional  night-time 
exc^tat^on  mechanisms  (i)  and  (ii) . 

In  addition  to  our  nominal  temperature  profile  Tu(z)  we  show  another 
temperature  model  TR(z)  on  fig. 5-2.  Thi<»  temperature  model  is  required  if  just 
the  mechanisms  (i)  and  (ii)  are  used  to  generate  a  calculation  that  will  fit 
the  data  D  on  Figure  5-1.  We  believe  it  is  very  unlikely  that  a  temperature 
profile  anything  at  all  like  T^(z)  actually  existed  on  4/11/74.  The  menopause 
temperature  data  provided  by  Ware  for  4/9/74  and  4/13/74  tend  to  substantiate 
our  contention. 

We  do  however,  believe  that  our  model  TN(z)  is  reasonably  realistic.  If 
this  is  the  case  then  there  must  be  some  CO^(\>^)  excitation  mechanism,  in 
addition  to  mechanisms  (i)  and  (ii),  which  is  operative  in  the  80  to  85  km 
region. 

We  considered  and  rejected  several  mechanisms  that  might  account  for  the 
additional  excitation.  •  For  example, energy  deposition  in  the  80  to  85  km  region 
by  hard  particle  precipitation  would  be  accompanied  by  visible  emissions. 

These  were  not  observed.  Or,  excitation  of  n|  via  the  reaction 

N  +  NO  -*•  N|  +  0 

would  require  unrealisticly  large  night  time  densities  of  N  in  order  to  provide 
the  required  source. 

The  only  mechanism  promising  enough  to  warrant  discussion  here  is  the  OH^ 
driven  mechanism,  that  is  production  of  N2  via  reaction  (5-1)  where  the  OH(v) 
is  initially  produced  via 


(5-3) 

(5-4) 

There  are  two  reasons  for  considering  this  mechanism:  (a)  The  85  km  peak 
altitude  requirement  could  be  nearly  satisfied  by  reaction  (5-1). 


H  +  03— 0H(v)  +  02 
or 

H02  +  0 — ►OH(v)  +  02 
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altitude  profile  for  0H  +  emission  typically  peaks  near  85  km,  (Baker  et  al., 

1977).  We  expect  then  that  the  altitude  where  peak  quenching  of  OH^  occurs 

2 

would  be  somewhat  lower  than  85  km.  (b)  The  total  power/cm  produced  by 

2 

reactions  (5-3) and  (5- 4) exceeds  the  power/cm  (apparently  of  the  order 
-9  2 

x  4  x  10  watts/cm  )  which  is  required  to  explain  the  "knee"  which 
occurs  near  85  km  in  the  11  April  1974  data. 

There  are  numerous  measurements  of  OH  airglow  volume  emission  rates 
available  in  the  literature  (Baker  et  al.,  1977).  To  investigate  the  OH(v)+N9  -*• 
n|  +  OH(v-l)  hypothesis  we  arbitrarily  selected  a  measured  airglow  OH  av  *  2 
volume  emission  rate  versus  altitude  e22(  v  "  2,z)  reported  in  the  literature 
by  Rogers  et  al. ,  (1973).  The  subscript  72  refers  to  the  date  these  data 
were  obtained  and  differentiates  them  from  a  2  cd  data  set  obtained  in  1975 
chat  we  discuss  later  in  this  paper.  The  volume  emission  rate 
published  by  Rogers  et  al.,  was  measured  on  6  March  1972  by  a  rocket-borne 
radiometer  with  a  passband  that  transmits  the  OH  Av  «  2  sequence  bands  4-*-2, 

3-+1,  and  2-+0.  To  estimate  the  corresponding  altitude  dependent  N*  production 
rate  QM*(z)  we  set  QN^lz)  *  e22(z)  <  k^/A  >  [n^J  *  The  calculation  <  k^/A  > 
requires  a  knowledge  of  the  rate  constants  for  reactions  (5-1)  and  (5-3).  For 
the  purpose  of  our  discussion  we  will  label  these  rate  constants  kp  k2>  and 
(kj  and  k2  refer  to  reaction  (5-1), they  are  defined  above).  A  knowledge  of 

V  ^ 

the  Einstein  coefficients  Av„  for  radiative  decay  of  OH(v')  to  OH(v")  is  also 
required.  For  our  calculation  we  used  the  A^„  calculated  by  Mies  (1974).  Mies' 
calculation  is  based  on  experimental  data  obtained  by  Murphy  (1971). 

In  the  initial  analysis  we  present  here  we  shall  neglect  production  of 
0H(v)  excitations  by  the  H02  +  0  reaction.  This  is  justified  since  we  know 
that  0^  +  H  -+■  0H(v)  +  0 2  contributes  significantly  to  production  of  the  night 
time  OH  airglow,  the  role  of  H02  +  0  in  contributing  to  OHfv)  excitation  is 
less  clearly  defined.  For  reaction (5-3) we  used  the  relative  0H(v)  formation  rates 
k3(9):  k3(8) :k3(7) :k3(6)  -  0.54:0.29:0.10:0.07  and  k3(v  *  5)  -  0  as  determined 
by  Good  (1976)  from  the  data  of  Charters  et  al. , (1976) . 

Quenching  rates  of  0H(v)  by  N2  and  02  for  4  *  v  «  9  are  given  by  Streit  and 
Johnson  (1976).  We  should  note  that  Streit  and  Johnson  were  not  able  to 
determine  if  the  generally  exothermic  reaction  5-1  proceeds  by  vibration  trans¬ 
fer  w  accompanied  by  partition  of  the  balance  of  the  energy  into  translation 
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and  rotation  vtr  or  whether  it  is  purely  vtr.  We  will  proceed  under  the  assump¬ 
tion  that  reaction  5-1  is  w  and  that  excess  energy  is  balanced  via  vtr.  We 
designate  the  rate  constants  for  w  excitationof  for  reaction  5-1  by  k^(v). 

For  O2  the  process  is  not  relevant  to  this  discussion,  we  designate  the  rates 
k2(v).  We  show  the  k^(v)  and  k2(v)  that  we  used  for  this  calculation  on  Fig.  5-3 
If  we  use  the  smallest  imaginable  values  k . (v)  -  0  for  i  «  1  and  2  and  for  v  ■  3 

ft  2  1 

we  get  E  e  I  dz  QN*  -  0.08  ergs/cm  sec,  and  if  we  use  the  largest  reasonably 

vv  J  *  -13  3,  < 

imaginable  values  k^(v)  -  10  cm  /sec  for  i  -  1  and  2  and  for  v  ■  3  we  get 

2 

-  0.19  ergs/cm  sec.  If  we  use  reasonable  Intermediate  extrapolations  as 

V  V  < 

shown  on  Fig. 5-3for  the  k^(v)  for  i  ■  1  and  2  and  for  v  ■  3  we  get  0.11 

ergs/cm2sec.  The  z  dependence  of  Qn|(z)  was  not  affected  appreciably  by  the 

extrapolation  assumptions.  This  point  is  illustrated  by  Fig.5-4where  the 

i  — 13  <  ^ 

profile  QN2(#  which  applies  for  k^(v)  ■  10  /sec  for  v  •  3,  the  profile  QN2 

which  applies  for  k .  (v)  for  v  -  3  extrapolated  as  shown  on  Fig.  5-3  and  the 
*  1  < 
profile  QN2  (-)  which  applies  for  kj (v)  «  0  for  v  »  3  are  presented. 

For  our  calculation,  we  shall  use 
QN*  which  we  designate  by  Qn|  72  in  the  remainder  of  our  paper.  This  will  be 
sufficient  for  our  purpose  of  demonstrating  that  the  mechanism  initiated  by 
reaction  (5-l)is  a  solid  candidate  for  explanation  of  the  enhanced  4.3  pm  zenith 
radiance  observed  on  4/11/74  in  the  80  to  85  km  altitude  region. 

For  the  reaction 


0H(v)  +  0  +  H  +  02  (5-5) 

-113 

we  use  k^  *  3.3  x  10  cm  / sec  (Westenberg  et  al.,  1970).  Although  k^  may 
have  a  small  v  dependence  (Rogers  et  al.,  1973)  this  is  not  particularly 
important  for  our  problem  since  reaction  (5-5)quenches  the  low  lying  0H(v) 
levels  most  efficiently  where  as  reaction (5-1)  produces  N2  mainly  from  the 
upper  levels  v  —  4  as  shown  by  Fig.  3-4. 

Other  quantities  necessary  for  the  calculation  of  <k^/A>  include  the 
0  and  N2  number  densities,  ^oj  and  respectively.  We  obtained  these  from 

the  DNA  Handbook  (Blank  et  al.,  1974).  The  volume  Emission  rate  e72  in  the 
bands  4+2,  3+1,  and  14-0  which  was  reported  by  Rogers  et  al,  is  also  shown  on 


■  ■■ — frjfcgg 
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Fig.  5-3.  The  data 
Johnston  (1976)  are  si 
Effects  of  variations 


X 


OH(V)+M  —►  OH(V-I) 


<8> 


r 


k 


Meos.  ex. 

•  ® 

X  ® 


JM _ 

02,  k2 
N2,  k| 


Measured  k  from  Streit  and 
Johnson,  J.Chem.Phys.64,95  1976 

i _ I _ I _ I _ I _ I _ i  i 

3  5  7  9 


or  k^v)  for  i  =  1  and  2  and  for  4$  v  <9  from  Streit  and 
own  here.  Reasonable  extrapolation  for  v  <  4  are  also  shown. 

in  the  extrapolated  values  are  discussed  in  the  text. 


ALTITUDE  (km) 


Fig.  5-4.  The  6  March  1972  measured  (Rogers  et  al)  altitude  profile  of 
volume  emission  in  the  bands  4  -*■  2,  3  -*■  1,  and  2  -*■  0,  and  the  inferred  vibra¬ 
tion  transfer  0H(v)  ^  N?  rate  QN?  are  shown  in  this  figure.  The  rates 
+  +  ±  *■ 

QNjM,  QN£  and  QN^-)  are  computed  on  three  different  assumptions  for  the 
value  of  the  rate  constants  k^(v)  for  v  <  3.  We  consider  the  rate  QN*  to  be 
the  most  reasonable  of  the  three. 


Fig.  5-4.  The  calculation  of  <k  /A>  which  is  appropriate  for  estimating  the 

*  ± 

transfer  rate  of  vibration  from  OH(v)  to  from  the  volume  emission  rate 
e22(z)  measured  in  bands  4-»2,-  3-*l  and  2-+0  by  Rogers  et  al.  (1973)  proceeds 
as  follows. 


<  VA ' £ 


N2]  “  nvkl(v>  W Vv- 


where  the 


nv  *  ^OH(v)j 


are  the  relative  populations  of  OH(v) .  Fojr  v  ■  9 


ny  ■  k3(v)  /  (A  +  Q^) 


(5-6) 


where  +  k$  [o]  ,  -  kl(v)[N2]  +  k2(v)  [o?]  and  where 

v  v  * 

A  ■  2-# A  »  •  For  successively  smaller  values  of  v  the  n  are  given  by 

« v  v 


nv  -  ^k3(v)  +  QB(v+1)Vl  +  V,)/-  QTv)  •  ^-r. 

Standard  C02(v3)  radiation  transport  methods  (Kumer  and  James,  1974;  Kumer , 
1977a)  may  be  used  to  compute  the  C02  4.3um  zenith  radiance  profile  R  that 
would  be  generated  by  a  source  of  given  by  QNf>7r  The  earths  atmosphere 
Is  optically  thick  for  C02  4.3um  in  the  80  to  loo’km  region.  For  this  reason 

there  is  a  maximum  in  the  zenith  radiance  profile  that  occurs  several  km  below 
the  maximum  in  Q  fff. 

The  altitude  of  the  maximum  in  the  OH  airglow  layer  is  known  to  vary  con¬ 
siderably  from  one  measurement  to  another  (Baker  et  al.,  1977).  It  would  be 
fortuitous  if  we  could  achieve  a  good  fit  to  the  data  (curve  D  on  Fig.  5-1) 
linear  combination  of  the  radiance  due  to  mechanism  (i)  and  (ii)  (Curve  A  on 
Fig.  5-1)  and  the  radiance  ROH72  computed  from  a  particular  set  of  OH  volume 
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emission  data,  such  as  the  data  «72(z)  rePorted  by  Rogers  et  al.  The  maximum  in 
the  Rogers  et  al.,  OH  volume  emission  data  occurs  at  a  relatively  low  altitude 
z  ■  84  km.  We  selected  a  2  cd  set  of  OH  volume  emission  data  c^(z)  with 
maximum  occurring  at  relatively  high  altitude  z  ■  89  km,  these  data  were 
obtained  3/4/75  (Baker  et  al.,  1977).  These  data  were  obtained  in  the  8-*8 
and  7-+5  bands.  The  vibration  transfer  rate  QN^  ^  to  is  given  by 


QN 


2,75 


'75 


I? 

V”  1 


\kl<v) 


K]/ 


v*7,8  nvAv-2 


and 

QN 


d  J dz  QN*  ■  0.0374  ergs / 
*  as  discussed  above. 


cm  sec.  The  radiance  Rqh7,.  “ay  bo  computed  from 


2,75 


A  least  squares  best  fit  to  the  data  D  of  radiance  curves  A,  Rqjj72  a°d 
ROH75  *8  shown  on  Fi8-  The  linear  combination  1.28A  +  0.816  RqH72  +  1*®* 
ROH75  *s  re<lu*re<*  to  achieve  the  fit  shown  on  Fig.  5-5.  The  coefficient  of 
curve  A  indicates  the  absolute  calibration  of  the  instrument  should  be 
modified  to  reduce  the  quoted  value  of  the  data  (  curve  D  on  Fig.  5-1)  by  a 
factor  1/1.28.  This  is  well  within  the  ±85%  quoted  accuracy  in  absolute 
calibration.  The  "recalibration"  is  justified  since  we  expect  our  calculation 
and  the  data  to  match  in  the  altitude  region  z  £  60  km  where  the  kinetic 
temperature  driven  mechanisms  (i)  and  (ii)  dominate,  and  where  we  have  ASL 
met.  rocket  sounding  data  on  atmospheric  temperature.  The  total  column 
vibrational  transfer  of  energy  from  OH*  to  N*  then  is  given  by 

-  (0.816  x  J  dz  QN*  ?2  +  1.84  x  I  dz  QN*  ?5> /l . 28 


*  .12  ergs/cm  sec. 


The  ±25°  uncertainty  in  mesopause  temperature  that  is  indicated  on  fig.  5-2 
causes  an  associated  contribution  to  uncertainty  in  the  recalibration  factor 
of  approximately  15%  and  to  uncertainty  in  of  approximately  ±17  Z. 

Also  we  should  note  that  the  9-*8  and  8-*  7  OH  bands  are  centered  near  4.47 

V  * 

and  4.14  \m>  Volume  emission  ev„  in  a  given  band  may  be  calculated  via 
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W/cm  sr 


Fig.  5-5.  A  least  squares  fit  to  the  data  D  which  is  given  by  1.28A  +  0.816 
+  1.84  R,,,,.,*  is  shown  on  this  figure. 

Ull 


85  89 

ALTITUDE  (km) 


Fig.  5-6.  The  composite  QN*  which  gives  the  fit  to  the  data  as  shown  on 
Figure  5-5.  The  composite  is  a  linear  combination  of  the  QN^  inferred  from  OH 
volume  emissions  antl  £75  that  were  measured  on  6  March  1972  and  4  March 
1975.  The  coefficients  are  0.816  and  1.84  respectively. 


V  =  I-816  e72  nv'Av(l/l!  VV 


v=2  v-2 


+  1.84  £  VAv"/  L  v' 

75  V  v-7,8  v 


Ay_2  ) / 1  * 2® 


The  calculations  shows  that  J"  dz  QN*  j  J  dz  j 


9  8 

e8  +  e7 


16.4  hence  direct  OH 


emission  near  4.3  pm  is  not  significant  compared  to  the  CO^v^)  emission  which 
results  from  vibration  transfer  from  OH(v)  to  N2<  The  fact  that  OH  v-*v  -1 
bands  which  lie  between  the  2.7  and  4.3  pm  features  are  not  observed  in  the 
CVF  spectra  supports  this  conclusion. 

+  vv 

On  Fig. 5-7we  show  how  the  OH^  ■+  N^  mechanism  affects  the  11  April  1974 
night  time  vibrational  temperatures  of  the  C02  001  state  T^^  and  of  N2  Tn2> 

The  solid  curve  is  nominal  atmospheric  kinetic  temperature  model  Tn>  The 

dash  dot  curve  shows  that  T^2  and  T^q^  are  identical  for  altitudes  z  <  ?5  km.  At 

z  >  72  km  the  dash  cot  curves  split  into  a  component  which  is  due  to  mechanisms 

(i)  and  (ii)  but  which  neglects  the  OH*  ^  N2  mechanism,  and  another 

i  w 

component  which  is  due  to  mechanisms  (i)  and  (ii)  and  the  OH  N2  mecha- 

t  vv 

nism.  It  is  apparent  that  the  OH  -►  N2  mechanism  dominates  the  T^and 
TQ01  in  the  region  z  >  73  km. 

We  show  Fig.5-7to  demonstrate  how  the  inclusion  of  the  OH*  ^  N^ 
mechanism  would  modify  the  calculational  results  for  night  time  vi!  rational 
temperatures  from  that  given  by  Kumer  and  James  (19/4).  The  corresponding 
modification  in  the  quiet  night  time  weak  C02  band  vibrational  temperatures 
may  be  calculated  by  methods  given  by  Kumer  (1977b) 

5.3  CONCLUSIONS 

We  conclude  chat  the  feature  near  85  km  which  is  observed  in  the  11  April 

%  w 

1974  4.3  pm  zenith  radiance  data  is  a  manifestation  of  the  OH  .*  N2  mecha- 
nism.  A  columnar  transfer  of  0.12  +  .025  ergs/cm  sec  of  vibration  energy  from 
OH*  to  N2  is  required  to  explain  the  feature.  A  carefully  planned  experiment 
utilizing  rocket-borne  measurements  of  OH  emission  and  C02  4.3  and  15  pm 
zenith  radiance  coordinated  with  simultaneous  and  independent  measurements  of 
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Fig.  5-7.  The  solid  curve  is  our  nominal  model  atmospheric  kinetic  tempera¬ 
ture  Tn  for  11  April  1974.  The  dash-dot  curves  are  the  coupled  N2  and  C02(001) 


vibrational  temperatures  TN  and  T^ .  The  dash-dot  curves  break  into  two 
components  for  z  -  72  km.  The  mechanisms  (i)  and  (ii)  were  used  in  the  compu¬ 


tation  of  both  the  dash-dot  curves  but  the  OH+  ¥*  N2  mechanism  was  included 


in  the  computation  of  the  upper  curve  but  not  in  the  computation  of  the  lower 


curve.  The  T^  and  Tftni  decouple  at  z  £  84  km  and  z  >  95  km  in  the  two  cases. 


001 


The  T^  are  the  dashed  curves,  the  Tq^  ^ottec*  curves.  This  figure 
illustrates  the  dominance  of  the  OH"  -*■  N"  mechanism  on  the  vibrational 


and  T„  above  z  =  72  km. 


temperatures  Tnm 


atmospheric  temperature  up  to  90  km  altitude  (by  ASL  sounding  rockets.  Nimbus 
6  PMR  atmospheric  temperature  sounder  and  ground-based  OR  rotational  temper¬ 
ature  measurements  for  example) ,  could  be  used  to  verify  the  role  of  the 
+  vv  1 

mechanism  OH-  N“  in  the  formation  of  4.3  pm  emission  in  the  altitude 
region  z  2  80  km. 
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6.0  ESTIMATES  OF  C02  HOT  BAND  EARTH LI MB  EMISSION  IN  THE  LWIR  REGION 

6.1  Introduction 

Sensitive  earthlimb  and  ATH  sensors  are  under  consideration.  Even  weak  back 
ground  emissions  might  impact  the  performance  of  such  sensors.  In  this 
chapter  we  report  estimates  fcr  the  CO^  LWIR  spectral  radiance  which  results 
from  auroral,  solar,  and  airglow  pumping  of  excited  states  of  CO^. 

6.2  Auroral  Processes 

The  00^(001)  will  be  produced  in  the  aurora  by  N*  ♦  CO^-^CO^tOOl) .  This  is 
followed  by  LWIR  emission  in  the  C02  laser  lines  via 

C02(001)  -C02(02°0)  + 

-CO2<10°0)  +hvwv. 

Subsequent  emission  will  also  occur  in  the  COj  LWIR  bands  2Uj+V2~*'2u2  +  ^2  " 


v'  v"  i' 
12  3 


U1  U2  V3 


2*1 


\  6m) 
Origins 


In  addition  LWIR  may  be  produced  by 

(l^O) 

CO. (001)  +  M-CO,  033  +  H  in  the  bands 

2  1  (oa1^ 


2y'  +  y'— 2vi  +  u2  "  3  —  2,  these  bands  are 


u'  u' 

v3  -  V1  U2  V" 

X  (um) 
Origins 

11*0 

0220 

13.48 

11*0 

02°0 

12.63 

11*0 

10°0 

14.52 

72 


V1  V2 

U3  ^ 

°1  i  °3 

a(m*) 

0330 

— 

0220 

Origins 

15.03 

0310 

— 

0220 

16.74 

0310 

- 

02°0 

15  45 

0310 

— 

10°0 

18.37 

6.3  Solar  Processes 

In  the  solar  case  the  CO^COOl)  Is  produced  by  absorption  of  4.3pm  sunllgnt  by 
C0_.  The  CO 2  LWIR  is  produced  from  the  C02(001)  as  discussed  above.  In  addi¬ 
tion  CO2(10°l)  and  CO2(02°l)  may  be  produced  by  absorption  of  2.7pm  sunlight. 
CO2  LWIR  is  generated  from  these  states  by  laser  radiation.. 


CO 


2 


02°l! 
10°1 
02  1 


12°0,04°0 

|12-0,04_0  +  hU(9  to  11pm). 

Il2  0,04*0 


The  lower  states  of  the  laser  transition  form  the  upper  state  for  LWIR  transi¬ 
tion  2V^  +  U*-~2V*'  +  v"  m  4—3,  Some  of  the  most  notable  of  these  transi¬ 
tions  are  listed  below 

-A. 

20°0  — 0310  11.56 

04°0  — 11*0  21.21  . 

Also,  the  mechanism  M  +  C02  J2]]  — M  +  C02  2"i+  *2*  "  5*  V3  “  °l  w111 

produce  LWIR  in  the  transition  2'  -*2^+1^ ■  5—4,  notable  of  these  are 

Origins 

11.66 


1330  —  0420 


6.4  Airglow  Processes 

In  the  night  time  airglow  COjCOOl)  is  produced  by  vibration  transfer  OH(u)  + 
-f7  ^2  +  CO2  -f7  COjCOOl),  LWIR  results  from  the  CC^COOl)  by  processes 

discussed  above. 


6.5  Summary  and  Results 

In  summary,  we  have  pointed  out  6  distinct  groups  of  CO2  LWIR  band  complexes  that 
may  be  pumped  by  the  aurora,  sunlight,  or  the  OH  airglow.  Below  we  shall 
designate  the  total  limb  radiance  in  a  given  band  complex  by  R^.  The  subscript 
n  designates  the  band  complex  as  follows: 


Subscript 


Process  and  Band  Corrplcx 


1  2  U1  +  U2  2  U1  +  m  2  **  1 
2  "  -  3  -  2 


3  ,r  4-3 

4  "  5-4 


5 


hu2.7  +  C°2 


101 

021 


+  laser  radiation,  9  to  11  nm 


6 


001 


020 

100 


II 


It  is  possible  to  use  the  altitude  profiles  of  vibrational  emperature  of  the 
CO2  upper  states  00°1,  02°1,  02^1,  and  10°1  in  order  to  estimate  (X^  LWIR  limb 
radiances  in  the  band  complexes  n  *  1  to  6.  In  order  to  calculate  the  limb 
radiance  estimates  we  used  CO2  vibrational  temperatures  Tv(z)  that  we  have 
calculated  for  the  March  24,  1973  auroral  breakup,  these  are  shown  in  Table  6-1. 
For  the  60°  solar  zenith  sunlit  case  we  use  the  Tv(z)  published  in  the  Kumer 
(1977b)  paper  and  for  the  airglow  case  we  use  Tqq^(z)  as  shown  on  Fig.  5-7  in  th 
the  previous  chapter.  The  limb  radiance  R  in  an  optically  thin  band  is 
approximately  given  by 

1  -  -T  /T  (z) 

=  Vn  S  a[C02]  e  0  7 


7  4 


R 


where  the  integral  is  along  a  limb  viewing  path  s,  A  is  the  Einstein  radiative 
decay  coefficient  of  the  upper  state,  [CO2]  the  CO2  number  density,  Tv(z)  the 
vibrational  temperature  in  the  upper  state,  and  To  =  he  E//k  where  E '  is  the 
vibrational  energy  of  the  tipper  state  in  units  cm  \  The  band  conplex  2  -•  1 
is  optically  thick  and  requires  a  more  conplicated  confutation 

1  „  -Tn /-  (v) 

Rq  =  ^  J  a  [C02]  ds  T  (aN(s))  Ae 

The  assumption  is  made  that  the  sub  levels  in  a  given  coup  lex 
of  levels  2  +  Uj  are  populated  according  to  local  thermodynamic  equilibrium 

(LTE).  This  is  probably  not  the  actual  case,  more  research  is  needed  to 
determine  the  population  distribution  within  a  given  conplex  of  levels  since 
this  determines  the  spectral  distribution  of  the  CO 2  LWIR  limb  radiance.  For 
present  purposes  we  will  use  the  estimates  based  on  the  LTE  assumption  in 
order  to  get  some  idea  of  the  C02  LWIR  limb  radiance  in  window  regions  that  is 
due  to  solar,  auroral,  and  airglow  punping.  The  estimated  radiance  in  each 
band  coup  lex  that  is  due  to  each  of  these  mechanisms  is  given  in  Table  6-2. 

It  is  possible  to  use  these  estimates  for  in  order  to  generate  synthetic 
LWIR  limb  spectra  for  tangent  altitudes  80,  90,  and  100  km.  Uiese  spectra  are 
generated  at  a  C Vi-like  resolution  of  .01  X  in  the  10  to  14.5  and  in  the 
16.5  to  20  nm  region,  and  spectra  with  0.5  cm  ^  HIRIS-like  resolution  are 
generated  in  the  region  of  the  Q  branches  of  selected  bands.  These  3pectra 
are  generated  by  line  by  line  methods  described  by  Kumer  in  HAES  Report  Mo.  70 
(1977).  Line  positions  and  strengths  from  the  McClatchey  tape  are 

utilized.  The  synthetic  3pectra  are  shown  in  Appendix  E. 
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7.0  EVIDENCE  FOR  A  PROMPT  4.3pm  AURORAL  EMISSION 


A  rocket  carrying  a  two  channel  liquid  nitrogen  cooled  radiometer  at  2.7  and 
4.3pm  and  a  photometer  at  4278  X  was  launched  10/26/78  at  the  Poker  Flat  Research 
Range,  Alaska,  for  the  Defense  Nuclear  Agency.  The  radiometers  sensitivity 
was  about  a  factor  X50  better  than  had  been  previously  used  for  any  rocket 
borne  observations  of  aurora  at  2. 7  and  4.3pm.  An  auroral  arc  (IBC  II+)  moved 
about  ten  minutes  before  launch  from  an  overhead  position  at  the  launch  area 
to  a  stable  position  about  50  km  to  the  north.  Comparison  of  the  4.3pm 
radiance  on  rocket  ascent  and  descent  at  altitudes  below  85  km  show  similar 
profiles  that  can  be  identified  with  the  OH  airglow  mechanism 

0H(v)  +  N0  V24  OH(v-l)  +  N* 

(L  4. 

followed  by 

4  *  C02  n  *2  -  -  c02  +  »Vj  m. 

From  90  to  about  105  km  the  ascent  zenith  radiances  are  nearly  an  order  of 
magnitude  brighter  than  on  rocket  descent.  The  ascent  90  to  105  km  altitude 
region  was  heavily  predosed  in  comparison  with  the  corresponding  region  on 
descent.  The  higher  ascent  values  of  4.3pm  zenith  radiance  are  consistent 
with  the  well  documented  (Kumer,  1977a),  long  time  response  4.3pm  auroral 
mechanism 

4  *  C02  **  K2  *  CO^COO!)  -  C0£  ♦  hX}  m 

where  the  ^  is  initially  produced  by  the  electron  precipitation. 


Ac  altitudes  above  about  100  km  for  both  ascent  and  descent  the  4.3pm  zenith 


radiance  data  show  structure  and  values  that  suggest  a  second,  nearly 
"prompt"  auroral  mechanism-  This  structure  is  clearly  evident  at  radiance 
levels  well  below  the  NER  of  previous  instrumentation. 

The  apparent  energy  efficiency  for  the  mechanism  is  t“0.3X  (or  1.1  x  10  ^ 
w/cm^  ster  at  A. 3pm  per  kR  at  A278  X)  for  altitude  Z  2  115  km.  The  data 
indicate  that  the  effective  value  of  c  decreases  rapidly  with  decreasing 
altitude  in  the  region  100<  Z  <115  km,  thus  suggesting  a  strong  quenching 
effect  in  this  altitude  region. 

A  preliminary  analysis  of  these  data  has  been  presented  at  the  DNA  IR  review 
meeting  held  10-12  April  1979  at  the  Air  Force  Geophysics  Laboratory  and  is 
summarized  in  the  document  AFGL-TM-18.  A  slightly  more  sophisticated  presen¬ 
tation  of  this  material  was  given  May  28  -  June  1  at  the  Annual  Spring 
Meeting  of  Americm  Geophysical  Union.  Vufoils  shown  at  this  presentation 
are  included  as  Appendix  F.  This  presentation  is  currently  being  written 
up  for  submission  for  publication  in  J.  Geophys.  Res.  An  important  goal  of 
future  field  experiments  ard  data  evaluation  is  of  course  to  attempt  to 
confirm  and  to  identify  the  mechanism  for  the  "prompt"  4.3um  aurora  so  that 
it  may  be  included  in  nuclear  IR  predictive  codes  such  as  ROSCOE.  Future 
rocket  borne  auroral  experiments  should  employ  high  spectra?  resolution  to 
identify  the  emitting  specie,  ard  should  be  designed  to  provide  information 
on  the  altitude  dependence  of  the  effective  4.3um  fluorescent  efficiency 
The  ground  based  support  should  be  designed  to  complement  the  rocket  data  in 
achieving  this  latter  goal  as  well  as  to  provide  input  for  a  reasonably 
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accurate  determination  of  the  time  response  of  the  nearly  prompt  mechanism. 

O  O 

Calibrated  ground  based  all  sky  red  (6300  A)  and  blue  (4278  A)  filtered  TV 
data  is  ideal  for  this  purpose. 

We  also  suggest  that  detailed  radiation  transport  analysis  in  application  to 
the  25  February  1974  side  looking  radiometer  data  may  reveal  confirming  evi¬ 
dence  for  the  existence  of  a  nearly  prompt  mechanism.  It  is  possible  that 
a  2  dimensional  analysis  of  the  24  February  1974  side  looling  data  might  also 
distinguish  the  contribution  of  the  dominant  slow  CO^  4.3um  aurora  from  that 
of  the  nearly  "prompt"  mechanism. 


8.0  CONCLUSIONS 


Methods  to  deal  with  the  time  dependent  two  dimensional  radiation  transport 
In  the  CO^  4.3pm  auroral  arc  were  developed  and  applied  to  a  nominal  aurora 
In  the  presence  of  a  wind  field.  The  nominal  auroral  arc  moved  In  from  the 
north  at  3  km/min,  and  then  stabilized  for  two  minutes.  The  calculations 
showed  that  lateral  gradients  in  the  nadir  radiance  that  are  predicted  if 
lateral  transport  is  neglected  are  considerably  steeper  than  for  the  predic¬ 
tions  which  do  include  the  effects  of  lateral  transport.  The  predicted 
smoothing  effect  of  lateral  radiation  transport  of  CO^  4.3pm  photons  in  the 
earths  optically  thick  upper  atmosphere  is  dramatic  enough  so  that  it  renders 
the  lateral  gradients  in  Che  CO^  4.3pm  auroral  nadir  radiance  no  more  severe 
than  the  lateral  gradients  in  auroral  2.7pm  nadir  radiance  even  though  the 
DC  level  of  the  C02  4.3pm  auroral  radiance  may  be  as  much  as  10  times  that 
of  the  2.7pm  aurora  in  steady  state  conditions. 

This  point  is  important  for  assessing  high  altitude  nuclear  degradation  of  the 
performance  of  advanced  nadir  viewing  scanning  sensors  due  to  high  altitude 
nuclear  mechanisms  such  as  x-ray,  (5-patch,  and  debris  patch  deposition  which 
are  similar  to  the  aurora.  To  provide  verification  we  considered  4  distinct 
experimental  concepts  applicable  for  the  4.3pm  region,  and  briefly  considered 
a  2.7pm  experimental  concept  too.  Most  of  these  concepts  involved  firing  a 
rocket  over  an  auroral  arc  and  viewing  the  aurora  in  the  nadir  in  order  to 
measure  lateral  structure  in  auroral  nadir  radiance  in  the  4.3  and  2.7pm 
spectral  regions.  Nadir  viewing  4.3pm  concepts  included  a  broad  b^and  filtered 
radiometer  (BBFR) ,  a  narrow  band  tilting  10  transmission  channel  picket  fence 
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filtered  iradiometer  (NBTFR)  anti  the  use  of  HIRIS.  In  each  case  the  instru¬ 


ments  are  to  be  used  with  a  wheel  of  C02  absorption  cells  in  order  to  provide 

-3  -1 

the  added  high  resolution  £6  x  10  cm  near  the  centers  of  the  CO2  lines 
that  is  necessary  to  discriminate  CO^  4.3pm  narrow  line  auroral  nadir 
radiance  from  the  earth  nadir  radiance  background.  Also,  0.2  to  0.5  cm  ^ 
resolution  between  the  lines  is  necessary  in  order  to  suppress  deleterious 
effects  of  the  broad  band  LTE  thermal  stratospheric  contribution  to  the  earth's 
4.3pm  nadir  radiance  background.  This  requirement  is  sufficient  to  rule 
out  the  BBFR  for  use  in  the  desired  arc  overfly  experiment.  The  NBTFR  and 
HIRIS  easily  satisfy  this  0.2  to  0.5  cm  1  spectral  requirement.  It  turns  out 
that  viewing  the  250°  earth  degrades  the  HIRIS  sensitivity  to  the  point  that 
verification  of  the  lateral  radiation  transport  smoothing  effect  on  the  C.O^ 
4.3pm  aurora  in  the  spatial  wing  of  the  nominal  arc  is  marginal.  However, 
since  the  NBTFR  views  the  250°K  earth  through  a  77 °K  filter  with  total  band 
pass  of  VI  cm  1  near  to  4.3pm,  sensitivity  is  not  degraded  appreciably  beyond 
the  inherent  InSb  detector  noise  limit  at  77°K  and  therefore  verification  by 
means  of  the  NBTFR  sensor  appears  to  be  eminently  feasible. 

The  optimum  design  NBFTR  uses  liquid  Nitrogen  (LN^,)  cooled  optics 

including  the  tilting  filters,  and  LN2  cooled  baffling.  This  reduces 

the  NEP  of  the  .71  cm  diameter  In  Sb  detector  to  =  3  x  10  ^  vh/ttz, 
nearly  the  inherent  LN2  cooled  detector  noise  limit.  The  narrow  band 

filter  is  a  4"  diameter  Silicon  (index  of  refraction  =3.4)  substrated 

Fabry  Perot  etalon  designed  so  that  the  centers  of  10  of  its  transmission 

1.2  16 

channels  roughly  coincide  with  the  lines  p8,  plO,  ...,  p26  of  the  C  02 
101  ♦  000  band.  A  300A  blocker  filter  allows  just  the  10  channels  to  be 
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transmitted.  The  narrow  band  filter  is  designed  so  that  the  transmission 

channel  full  width  half  maximum  FWHM  “0.21  cm  *  (i.e.,  4A)  for  an  FOV  with 

the  1/2  cone  angle  a  -  2°.  The  through  put  of  the  optimum  NBTFR  is  “.062  cm2 

-13  2 

sr  for  an  overall  optics  transmission  of  “20Z.  The  NER  “5.0  x  10  w/cm  sr. 

Spectral  scanning  between  the  CO2  line  centers  is  achieved  by  tilting  the 

‘  -3  -1 

filter.  High  spectral  resolution  of  “6  x  10  cm  near  the  CO2  line  center 
is  achieved  by  utilization  of  the  CC^  absorption  cells. 

As  part  of  the  study  it  was  also  verified  that  an  optimum  NBTFR  could  be  used 
effectively  in  an  arc  overfly  experiment  to  verify  that  lateral  gradients  in 
the  2.7pm  auroral  nadir  radiance  are  indeed  as  sharp  and  dramatic  as  one  would 
expect  from  an  essentially  prompt,  optically  thin  auroral  emission. 

Estimates  for  the  spatial  structure  to  be  observed  by  the  sensor  module  on  the 
spatial  EXCEDE  daughter  vehicle  on  viewing  the  electron  gun  discharge  at  the 
mother  vehicle  at  90  km  altitude  were  also  generated.  These  suggest  that 
prompt  emission  by  NO+  might  be  expected  to  dominate  the  region  of  instantaneous 
energy  deposition  but  that  the  long  lived  mechanism  N^"  +  CO^-*-  CO^COOl)-* 

+  hv  4.3pm  will  dominate  the  predosed  regions.  Also  it  is  estimated  that 
radiation  transport  into  regions  that  are  adjacent  to  the  region  of  instanta¬ 
neous  deposit  and  to  the  predosed  region  will  be  strong  enough  to  produce  a 
signal  above  background.  Data  that  will  be  obtained  in  these  so  called  transport 
regions  can  serve  as  a  guide  for  the  development  of  a  numerically  effective 
CO^  4.3pm  3  dimensional  radiation  transport  facility  for  application  to  the 
nuclear  case  (highly  structured  8  deposition  region). 


Analysis  of  a  feature  near  85  km  which  is  observed  in  the  11  April  1974  4.3ptn 

zenith  radiance  data  indicates  that  it  is  a  manifestation  of  OH  airglow  related 

mechanism  OH  for  producing  4.3pm  and  LWIR  emission.  A  columnar 

2 

transfer  of  0.12  +  .025  ergs/cm  sec.  of  vibration  energy  from  OH  to  is 

required  to  explain  the  feature.  Although,  to  date,  there  has  been  no  experi¬ 
ment  directly  aimed  at  simultaneous  measurement  of  all  the  parameters 
(temperature,  OH  emission,  4.3pm  and  LWIR  emissions)  that  are  pertinent 
to  the  OH  airglow  mechanism  for  CO^  4.3pm  and  LWIR  emission,  still  there  is 
considerable  evidence  (Kumer,  1976;  Kumer  et  al  1978,  this  report,  etc.)  for 
the  existence  of  the  mechanism.  Detailed  verification  and  quantization  of  the 
mechanism  could  be  achieved  for  examply  by  carefully  planned  experiment 
utilizing  rocket-borne  measurements  of  OH  emission  and  CO^  4.3  and  15pm  zenith 
radiance  coordinated  with  simultaneous  and  independent  measurements  of 
atmospheric  temperature  up  to  90  km  altitude  (by  ASL  sounding  rockets.  Nimbus 
6  PMR  atmospheric  temperature  sounder  and  ground-based  OH  rotational  tempera¬ 
ture  measurements  for  example). 

The  mechanisms  aurora,  airglow,  and  sunlight  fluorescence,  for  pumping  CO^ 
emissions  in  LWIR  earthlimb  window  were  addressed.  These  mechanisms  are  impor¬ 
tant  in  the  natural  background  and  have  nuclear  analogues.  Estimated  synthetic 
spectra  produced  by  these  mechanisms  were  generated  for  tangent  altitudes  80, 
90,  and  100  km  at  a  CVF-like  resolution  of  .01  A  in  the  10  to  14.5um  region 
and  In  the  16.5  to  20pm  region,  and  also  with  0.5  cm  *  HIRlS-like  resolution 
in  the  region  of  the  Q  branches  of  selected  bands. 


In  many  cases  the  estimated  spectral  radiances  were  well  above  realizeable 
sensor  NESRs. 

A  major  uncertainty  in  these  estimates  is  that  they  were  calculated  on  the 
assumption  that  the  sub  levels  in  a  given  complex  of  CO^  levels  2  +  i>2  are 

populated  according  to  local  thermodynamic  equilibrium  (LTE) .  This  is  probably 
not  the  actual  case,  more  research  is  needed  to  determine  the  spectral  distri¬ 
bution  within  a  given  complex  of  levels  under  atmospheric  conditions  near 
100  km  altitude  since  this  determines  the  spectral  distribution  of  the  CC^ 

LWIR  limb  radiance. 

The  analysis  of  4.3pm  auroral  zenith  radiance  data  obtained  10/26/78  via  a 
rocket  borne  radiometer  of  improved  sensitivity  indicated,  for  the  first 
time,  the  presence  of  a  relatively  weak  mechanism  for  prompt  auroral  emission 

at  4.3pm.  The  apparent  energy  efficiency  for  the  mechanism  is  e~0.3%  (or 

-10  2  ° 

1.1  x  10  w/cm  ster  at  4.3pm  per  kR  at  4278  A)  for  altitude  Z  £  115  km. 

The  data  indicate  that  the  effective  value  of  e  decreases  rapidly  with  de¬ 
creasing  altitude  in  the  region  100  <  Z  <  115  km,  thus  suggesting  a  strong 
quenching  effect  in  this  altitude  region. 

An  important  goal  of  future  field  experiments  and  data  ecaluation  is  to  attempt 
to  confirm  and  to  identify  the  mechanism  for  the  "prompt"  4.3pm  aurora  so  that 
it  may  be  included  in  nuclear  IR  predictive  codes  such  as  ROSCOE.  Future 
rocket  borne  auroral  experiments  should  employ  high  spectral  resolution  to 
identify  the  emitting  specie,  and  should  be  designed  to  provide  information 
on  the  altitude  dependence  of  th  effective  4.3pm  fluorescent  efficiency. 

The  ground  based  support  should  be  designed  to  complement  the  rocket  data  in 
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achieving  this  latter  goal  as  well  as  to  provide  input  for  a  reasonably 
accurate  determination  of  the  time  response  of  the  nearly  prompt  mechanism. 

O  O 

Calibrated  ground  based  all  sky  red  (6300  A)  and  blue  (4278  A)  filtered  TV 
data  is  ideal  for  this  purpose. 

We  also  suggest  that  detailed  radiation  transport  analysis  in  application  to 
the  25  February  1974  side  looking  radiometer  data  may  reveal  confirming  evi¬ 
dence  for  the  existence  of  a  nearly  prompt  mechanism.  It  is  possible  that 
a  2  dimensional  analysis  of  the  25  February  1974  side  looking  data  might 
also  distinguish  the  contribution  of  the  dominant  slow  CC^  4.3pm  aurora  from 
that  of  the  nearly  "prompt"  mechanism. 
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Appendix  A:  The  Approximate  3D  Formulation 

Here  we  give  the  3 -dimensional  formulation  for  C02  4.3  uni  excitation  in  the 
earthfe  upper  atmosphere.  The  application  for  other  problems  involves  an  appro¬ 
priate  change  In  the  form  of  the  Green's  function  H^.  The  formulation  is 

S5D(a»,x,y)  -  SQ(aN,x,y) 

+  o  J  dN'dx'dy'fi(lOS3D  (or,x',y')^D(|B-r|,|x-x'|,|y-y'|)  ,  (al) 

where 

H3D-Hbfa'K)/lw(|?  -?'|)2  •  (a2) 

Approximate  forms  are 
SID(aN,x,y)  -  So(oN,x,y) 

+  a  J  dN'n(»')S1D(aN%x/y)J*  dx'dy'H^lN-N'f ,  |  x-x'| ,  |y-y'| )  (a3) 

and 

S^D(oN,x,y)  «  S^(oN,x,y) 

-♦  aJdx'dy^NjS^pCoNjx'^')!  dN'H5D(|N-N#|,|x-x'|,|y-y'j)  ,  (a4) 

where  . 

-  s1D(i  -  n(i-E)  ). 

Useful  identities  are 

I  dxdy  S3D  -  J*  S^D  “  I  ^  S1D  (a5) 

and 

°J  «  S5D  ‘  °I  “  ^D-  (“6) 
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Appendix  B:  Nominal  Arc  Model 


The  nominal  arc  was  assumed  to  move  out  of  the  north  »t  •  3  Ir  t/min.  and  it 

was  assumed  to  have  stopped  and  stabilized  its  motion  for  two  minutes  prior  to 

the  time  for  which  the  calculations  on  Figure  2-1  apply.  The  nominal  arc  spatial 

structure  is  modeled  on  the  basis  of  structure  observed  In  the  3/27/73  arc 

in  the  ICECAP  measurements  program.  In  Fig.  B1  we  show  the  deposition  at 

arc  center  qt,(z)  and  at  tie  point  of  rocket  penetration  q  (Z)  for  the 
H  up 

March  27,  1973  BB73  rocket  measurement  of  4.3pm  auroral  zenith  radiance.  The 
brightness  of  the  3914  X  aurora  is  95  and  15  kR  respectively  for  these  two 
deposition  profiles.  The  deposition  shapes  q^  and  q^p  are  discussed  in  more 
detail  by  Kumer  in  AFCRL-TR- 74-03 34.  In  order  to  form  a  nominal  aurora  of 
30  kR  3914  X  brightness,  105  km  altitude  peak  in  deposition,  and  15  km  width, 
we  utilize 


q(z,  0)  «  qH(z  ♦  5) 

and 

q(z,  xL)  -  —  qup  (z  -  5) 

where  x^  -  7  km  x  /in  (95/15)'  «  9.5  km.  The  deposition  then  in  the  arc  at 
altitude  z  and  position  x  is  given  by 

2 

/  \  /  -(x/o(z)) 

q(z.  x)  -  q(z,  0) e 


where  o(z)  -  X^/Zln^z) /qup(z) ) . 

Arc  motion  and  winds  are  accounted  for  as  follows 
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Figure  Bl.  Models  for  q(z)  at  the  upleg  penetration  point  q^Cz)  and  at  0X0  center  q^z)  for  the 
March  27,  1973  B373  rocket  shot  (see  AFCRL-TR-74-033^  for  more  discussion  of  a  and 


,  ,  n,  -(x’/a(z))2 

q(z,  x  )  =  q(z,  o)e 

whe  re 

x'  =  X-(VA  -  Vw(z)(t+2) 
for  t  <  -2  min.  before  stabilization  time  and 

i  ( ? ) 

x  =  x  +  '(t+2)  after  arc  stabilization  time.  The  wind  fields  V^z 

measured  by  tracking  a  TMA  release  on  the  3/23/73  date  by  Faire  and  Bettinger 
(private  communication,  1976)  are  shown  on  Figure  B2. 

The  So(oN,X)  that  are  shown  on  Figure  2-1  in  the  text  arc  calculated  by  the 
STDA  method  (Kumer,  Equation  (7)  in  AFCRL-TR-74-0334,  1974) 

o  ^  t  . 

S  (oN.x)  =  ef (z)  /  t  dt  et/T  q(z,x') 

° 

where  x(z)  is  the  CO^  4.3um  relaxation  time  given  by  the  curve  EFA,  Qw 
that  is  shown  in  Figure  1  of  the  publication  theory  of  the  CC^  4 . 3um  aurora 
and"Related  Phenomena"  by  Kumer  (1977a)  and  where  t  and  f(z)  are  defined  in 
the  report  AFCRL-TR-74-0334.  Here  we  used  t  =  10  Nn  vibration  quanta/ion  pair. 
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Figure  B2.  Preliminary  wind  speed  data  Kindly  provided  by  Ffcire,  Bettinger,  et  al. 


FT-T 


Appendix  C:  Procedure  I'cr  Rapid  Converge) ice 


Kquationa  such  as  eq.(2-l)  wuso  ic  discretized  to  a  matrix  formulation  (see  KJ) 
it*  order  to  obtain  a  solution.  The  matrix  formulation  of  eq.  (2-1) 

S  =  £  +  HS  (Cl) 

o 

or 

(l-H)S  =  SQ.  (C2) 


The  solution  is  then 


where 


S  =  GS 

o 


G  =  (1-H)'1. 


(C3) 


(C4) 


If  there  are  too  many  mesh  pointy  it  is  not  practical  to  solve  for  G.  In  this 

case,  we  must  use  an  iterative  technique.  We  use  a  trial  solution 

S'  a  S  +  6  where  S  and  6  are  unknown.  After  n  iterations  of  eq.  (Cl), 
we  get 

S  =  S  +  H°  6-  (C5) 

n 

For  any  physically  realistic  radiation  transport  H, 

h"  -»  0  as  n  -*  tc. 

Thus,  satisfactory  convergence  is  achieved  if  we  iterate  enough  times. 


If  H*1  -*  0  quite  slowly  as  n  increases,  it  could  be  expedient  to  attempt  an 
iteration  scheme  which  provides  more  rapid  convergence  with  increasing  n. 

For  example,  we  may  use  an  approximate  inverse  matrix  G'  given  by 

G'  =  G  -  K  .  «:?) 


where  G  and  K  are  unknown,  and  go  through  the  itera: i  n  scheme 
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UNCLASSIFIED 


S  -  S  +  HS'  -  S  ♦  iT5, 

F  o 

(l-H)fi  -  S'  -  Sp> 

Sx  -  S'  -  G'(l  -  H)6 
-  S  ♦  K6, 

and  ther\  after  n  iterations, 

S  =  S  +  K°  6. 
n 

If  K°  -*  0  appreciably  faster  than  H°  -*  0  with  increasing  n,  then  more  rapid 
convergence  is  achieved  if  we  use  the  iteration  scheme  given  by  eq. (C7) 
through  (cil)  rather  than  simple  iteration  of  eq.  (cl) . 

For  the  solution  of  equation  (2-1),  we  found  the  approximation 

G'  «=  1/(1-Hii) 

to  be  effective  for  achieving  rapid  convergence  via  the  interaction  scheme  of  eqns. 
(C7)  through  (Cll) .  With  this  approximation  for  G*  30  iterations  of  eq. (C7) 
through  (Cll)  were  required  to  achieve  satisfactory  convergence.  However,  130 
simple  iterations  of  eq.  (Cl)  were  required  to  attain  conyarable  con¬ 
vergence.  The  trial  solution  S°  *  was  used  in  both  cases.  The  scheme  of  eqns. 

(C7)  through  (cll)  saved  in  computer  time  a  factor  0.23  over  single 
iteration  of  eq. (Cl) . 


(C8) 

(C9) 

(CIO) 

(Cll) 
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In  order  to  achieve  satisfactory  convergence,  we  required  that 

,6  =  |  S  ,  -  S  I  S  10‘5  x  £ 

n  1  n+1  n<  o  max 

holds  for  all  elements  6  where  S  is  the  largest  element  of  the 

n  o  max 

inhomogeneous  vector  SQ.  This  is  reasonable  when  one  considers  that 

fdx6n  =  J  dx(K°  or  1^)6  *  0  is  an  identity  if  the  initial  trial  solution  is 

given  by  .  This  is  true  since  6  *  -  S  so  that  fdx  6  «  0  [by  the 

identity (2-5) in  the  text].  The  quantities  H*1  and  K11  must  be  functions  of 

|x  -  x'|  .  Next,  we  consider  that  Jdx  dx'g(|x  -  x'|)6(x')  »  Gjdx'  6(x')  ■  0  where 

G  ■  Jdx  g(  |  x  -  x'|  )  and,  therefore,  r  G  «  0.  Hence  Jdx  Hn6  *  J*dx  1^6  ■  0  » 

Jdx^n  -  S). 


APPENDIX  D 


CALCULATIONS  FOR  SPATIAL  EXCEDE 


Introduction: 
estimates  R*, 


In  this  section  we  provide  the  details  of  our  calculation  of 

_II  ,  _  that  are  shown  in  table  4-1  of  the  text. 

R  and  RN(>f 


Problem  Definition  for  CO^  Emission:  The  radiation  transport  problem  is 
essentially  3  dimensional  (3D).  As  yet  we  have  not  developed  calculational 
facility  to  deal  with  3D  problem,  however,  for  this  study  we  used  a  version 
of  the  3D  approximation  as  described  in  our  paper  "Approximate  and  Exact 
Technique  for  Multidimensional  Radiation  Transport  in  a  Plane  Parallel  Atmos¬ 
phere;  Application  to  the  4.3  ym  Auroral  Arc"  (Kumer,  J.  Ouant.  Spectroscp.  and 
Radiat.  Transfer  19,  649,  1978). 

We  approximated  the  deposition  by  the  rectangular  geometry  described  in  chapter  4. 
Next  we  generated  time  dependent  ID  solutions  S^R  (see  Appendix  A)  for  initial 
electron  gun  induced  vibrational  excitation  of  CO^  (call  it  sj^)  and  for  initial 
excitation  of  N^  (call  it  S**)  by  methods  described  in  our  paper  "Theory  of  the  CC>2 
4.3  ym  Aurora  and  Related  Phenomena"  (Kumer,  J.  Geophy.  Res.  82,  2203,  1977a). 

The  inhomogeneous  terms  which  are  required  to  generate  (these  are  defined  in 
the  paper  just  cited)  were 
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120  kw  x  (l  cv/l.6  E-22  kw)(lO  quanta/ionization) 

^1  "  (7E5  cm  x  2  E3  cm  x  1  cm)(35  ev /ionization) [R? ] 

on  the  interval  90  <  z  <  90.0  km,  and  n1  =  0  elsewhere,  *  0  everywhere. 
The  time  dependence  of  n1(x,z)  is  given  by  r^(z)  /  0  for  -  10  m  <  vt  <  10  m 
where  x  measures  horizontal  distance  perpendicular  to  the  magnetic  field  and 
v  =  0.6  km/sec.  In  order  to  compute  it  is  merely  necessary  to  interchange 
the  identities  of  and  r^.  Thus  we  are  assuming  that  the  electron  gun 

and  C02(v^)  are  identical  (i.e.,  «  ecQ  ). 

This  is  legitimate  for  the  purpose  of  estimating  the  EXCEDE  signature  for 
these  2  mechanisms.  Performance  of  the  experiment  should  ideally  result  in 
the  actual  determination  of  the  absolute  excitation  efficiencies  for  both 
species. 


excitation  efficiency  of 
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Some  representative  solutions  and  are  shown  in  table  D.l.  These 
solutions  show  that  forms  a  DC  trail  in  the  predosed  region  (FR)  Just  as 
expected,  and  that  emission  due  to  is  essentially  prompt,  again  as  is 
expected. 

In  addition  in  the  predosed  region  radiation  transport  produces  a  signal  in 
the  70  m  altitude  region  which  is  directly  below  the  70  m  predosed  region. 

This  transport  signal  is  reduced  by  about  a  factor  1C  from  that  in  the  pre¬ 
dosed  region  itself.  We  will  show  below  that  the  transport  signal  is  com¬ 
parable  with  the  ambient  radiance  level  one  expects  at  90  km.  In  addition 
there  is  a  noticable  gradient  across  this  JO  m  transport  region. 

Three  Dimensional  Transport  Effects:  Since  emission  in  the  dosed  and  predosed 
regions  is  essentially  the  first  emission  for  a  given  quantum,  the  effects  of 
3D  multiple  absorption  and  emission  (scattering)  in  the  dosed  and  predosed 
regions  will  be  negligible.  These  effects  will  be  substantial  in  the  transport 
region,  however,  since  emission  in  this  region  results  from  one  or  more 
scatterings.  Just  vertical  radiation  transport  is  accounted  for  in  the  cal¬ 
culation  of  however,  it  is  shown  in  appendix  A  that 

Jd  x  dy(S1D  -  S5D)  -  0 

where  S  is  the  exact  3D  radiation  transport  solution.  Radiance  values 
3D 

R(z,x)  viewing  from  the  side  are  computed  from  the  source  function  S(z,x,y)  via  an 
integral  along  the  y  direction  (1  to  both  the  magnetic  field  and  the 
velocity) 
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Table  D.  la  .One  dimensional  solution  S*  (z,x)  for  mechanism  I. 
Radiance  R*  In  units  watts/cn^sr  glvefrby  %  R*  —  o  [COblS  . 
Also,  we  note  here  that  In  order  to  perform  the  calculation  11 
was  numerically  convenient  to  use  z  ■  89*97  km  rather  than 
z  -  90  km  as  the  lower  boundary  of  the  deposition  region. 
However,  throughout  the  rest  of  the  text  z  ■  90  km  Is  taken 
as  the  lower  boundary. 


z(kml 


-  26  m,  -  1.5  km. 


167  km 


140,00  •1945*01  L J404-0I  ,1970-01^. 1569-00  .1054-09 

107,40  ,2440-01  ',4079-05  ,2641-01  ,4043-00  ,5044-09 

90,00  .3020-05  , 0114-05 _ ,4016-05  .1101-01  >0701-09 

95,12  .0454-05  1.1009-04  ,0404-05! ,5497-00  ,1304-00 

92,57  ,1405-04  !, 1905-04  ,7041-05  ,0149-00  ,1944-00 


rJ  k 1  Ti»r 


iinT 


mw 


lit  TT*T 


90,23  ,7490-04  ,0401-04  ,0455-05  .1914-07  .2763-06 
90,04  ,0015-03  ,0907-03  ,0297-05  ,5075-07  ,3944-00 


09*99  *9041^of  ,9199-03  1  ,0294-05  ,0173^07  ^3902-04 

09,95  ,1246-03  ,1360-03  ,0170-05  ,2500-07  ,2900-00 


liiniLi 


►xf^rr»w 


hil'in*! 


1000-03  ,1137-03  ,6104-051,2304-07  ,4919-00 
1016-04  ,1473-04  ,3333-051,1001-07  ,2503-00 
6396aQ6  ^.1033- Q5  ,5770-06  ,2409-00  ^147A-Q0 
9932-07  ,1312-06  ,6009-07  ,5042-09  ,6043-09 
2913-07  ,4203-07  , 1573-07| ,2149-09  ,3466-09 


niTri'ii 


mw  JI  T.TJ'l 


75.20  ,0964-00  ,1241-07  ’  ,4074-001 .8539-10  .1759-09 
71,34  ,1333-00  ,2039-00  ,3940-09! ,2199-10  ,4292-10 
44-9A — ^7-li3-40-i  .9415- IQ  ,2373-10  .2653-n  .3731-1 1 

37.21  .9267-11  ,1197-10  ,3346-11  ,6905-12  ,5126-12 
33,63  ,2564-11  ,3320-11  ,1065-111,2967-12  ,1350-12 


I  \  was 


50,26  ,7916-14  ,1030-11  ,3606-12  ,1372-12  .3S92-13 
50,00  .7209-14  ,9304-12  ,3364-12  ,1031-12  ,6000-14 
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Table  D.  16.  One  dimensional  solution  (*,x)  for  ^mechanism  II. 

Radiance  R**  is  given  sis  units  watts/cm  sr  by  %R  *  o[COg]S  . 


160.00  1 
107.10 

99.12 

92.97 


hrJ 


90.23 

90.04 

90.02 

69.99 

89.99 


X  »  0, 

•4746-00 

•9797-00 


-6  m,  -  26  n, 

•  1337-09.  •1913-09 
•1070-09  •7116-09 


■1.5  km,  -  167 


•1937-09 
,3229-09  I 


HiniLI 


•  2421-04  , 
.4201-03 


.4322-03 

.4390-04 


.4413-09 

.8909-05 


nuini 


,9071-04| 

•7147-03 

■7296-03 

.7309-03 

.8088-04 


,1616-04; 
. 2900-041 


■  in 


•1093-03 
•1079-02 
.1096-02 
•1111-02 
9 1700-03 


•0709-09 
•6807— Of 
♦1441-00 
•2204-04 
.3364-04 


•1115-03 

.1046-02 

,1009-02 

.1061-02 

•1711*03 


•1479-09 

•1923-09 

•1905-09 

.9177-09 

.7931-09 

•  ♦  i  1.  Am 


•1331-04 
•4091-04 
•  4009-04- 
•4092-04 
.1990-04 


50.20  .2997-12  .1333-11 

90,00  1  •2400-12  *9230-12  .1392. U 


•  1372-10  '  .1120*09 

•  1194-10  :  .2136-10 


%  R5D  -  fa  [C02]dy'n(z)T(otC02]|y-y#|)S3D 


so  that  in  the  case  of  S?f  where  ~  S._  ~  0  for  x  <  -  .01  km  we  therefore, 

ID  dx  ID- 

expect  that  R^^  R^  for  x  <  -  .01  km.  On  the  other  hand  since  is 

concentrated  in  the  instantaneous  deposition  region  (IDR)  we  expect  that  R^ 

for  -  .01  <  x  <  .01  km  and  that  R*  <  R*  on  the  rest  of  the  x  axis. 

ID  3D 


In  order  to  quantitatively  estimate  these  effects  we  attempted  to  calculate 

the  approximate  3D  solution  which  is  defined  in  Appendix  A.  Unfortunately 

the  scope  of  the  effort  did  not  permit  the  development  of  a  working  computer  program 

A 

for  the  calculation  of  S^D  for  application  to  EXCEDE  planning. 


We  did  ,  however,  develop  a  working  computer  program  which  uses  a  diffusive 
approach  to  the  problem.  This  approach  is  less  accurate  but  easier  to 
implement  than  is  calculation  of  by  solution  of  the  lateral  ?D  radiative 
transport  integral  equation  (aU)  in  appendix  A.  In  order  to  calculate 
the  diffusion  approximation  solution,  we  perform  the  following  calculation 


nA  D 
3D 


(x,y)  =  J  dx'J  dy'GD(f;,t)S1D(x')6(y#) 
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x  a  .01  km, 

o 


G 


D 


1  e-$2A*t 

%  at  ' 


2 


,  ,.2  2 

(x-x  )  +  y 


t  =  (xQ  -  x')/v. 


p 

The  value  a  *  .1*1  km  /sec  is  the  radiative  diffusion  coefficient  appropriate 
for  90  km  altitude.  The  derivation  of  this  value  is  given  in  the  report 
(Kumer,  J.  B. ,  "Theoretical  Estimates  for  the  Lateral 

Structure  in  the  4.3  um  Auroral  Arc,"  Draft  Report  In  Utah  State  University 
Subcontract  SC-75-019,  Authorized  by  Government  Contract  F19628-73-C-0048; 
February,  1976) .  This  approximation  conserves  the  identity 

Jdxdy  S5D  =  Jdxdy  S1D  «  Jdxdy  S^. 

Resurts  for  CO^  and  N,,  -»  CO,,  Radiances:  A  plot  of  the  radiance  component  RI(x) 
which  results  from  direct  electron  beam  excitation  of  COp(v^)  (mechanism  i) 
is  shown  on  fig.  D-l.  The  approximate  radiance  calculations  R*^  afld  R^  are 
shown.  These  are  computed  from  the  approximate  solutions  and  S^ 

respectively.  Radiance  values  for  z  ■  20  m  above  and  below  z  -  90  km  (which  is 
the  lower  boundary  of  deposition) are  shown  on  fig.  D-l. 


The  estimates  R^°  are  shown  by  solid  curves  and  R^  by  dashed  curves, 
label  A  designates  the  case  where  z-90km»+20m  and  B  the  case  where 
z  -  90  km  ■  -  20  m. 


The 
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Figure  D-l  The  Solid  Curves  are  our  Estimates  for  the  Mechanism  I 
Component  of  Lateral  4.3pm  the  Radiance  Profiles  R|  3D 
Taken  20  m  Above  (Case  A)  and  20  m  Below  (Case  B)  the 
Boundary  of  the  Dosed  and  Predosed  Regions.  The  dashed 
curve"  are  lateral  radiance  profiles  R^D  for  the  Same 
Circumstances.  These  curves  apply  for  deposition  at 
90  km.  The  difference  in  Cases  A  and  B  in  the  transport 
regions  (which  are  x  >  10  m  and  x  <  -  10  m)  is  an  arti¬ 
fice  of  our  approximate  solution  as  is  explained 

in  the  text.  We  expect  R^  —  10"^  watts  is  probably 
fairly  close  to  the  actual  value  which  should  obtain  in 
the  transport  regions. 
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x pec  ted  on  the  bar  is  of 


We  note  that  >  R*^D  in  the  region  IDR,  just  ns  e 
ur  earLier  discussion.  In  the  PDR  we  see  that  transport  dominates  over  pre¬ 
dosing  in  generating  R^.  Next  we  note  a  large  disparity  in  the  value  of  R^5 

in  the  region  ATR  and  PTR  as  compared  to  that  in  the  region  LTR.  Intuitively 

,IAD 


we  know  this  disparity  exists  in  the  calculations  because  use  of  the  S 


3D 


approximation  does  not  allow  for  a  direct  contribution  via  the  Integral  in 
IAD 

equation  2-1  to  in  the  region  LTR  from  the  initial  emission  which  occurs 

in  the  region  IDR.  However,  since  we  know  JdN(S^  -  S^)  =0  is  an  identity, 

it  seems  reasonable  that  we  take  R^  *  (?  x  10  ^  ;  x  10  10)/?  watts/cm^  sr  H 

-Q  ° 

10  watts/cm^  sr  to  obtain  near  to  the  region  IDR  in  all  the  transport 
regions,  PTR,  ATR  and  LTR. 


The  OH*  +  N0  -*  N* ;  N*  +  CO,^  -»  C0,^(v  )  generated  ambient  1.3  pm  radiance 
Level  for  horizontal  viewing  at  90  km,  as  inferred  from  the  1.3  pm  zenith 
radiance  data  set  which  was  obtained  on  l/ll/7l,  is  also  shown  on  fig.  D-l. 
l/ll/7l  ambient  90  km  horizontal  1.3  pm  radiance  thus  shown  is  probably  some¬ 
what  smaller  than  would  be  expected  on  the  average. 

It  could  conceivably  be  greater  by  a  factor  0  as  is  indicated  on  fig.  D-l.  How¬ 
ever,  we  can  in  principle  launch  the  EXCEDE  experiment  on  an  evening  for  which 
the  ground  based  OH  airglow  is  observed  to  be  lower  than  average,  and  hopefully 
in  this  way,  seLect  an  evening  for  which  the  ambient  -4.3  pm  radiance  is  lower 
than  usual,  as  was  the  case  on  the  evening  of  ■«/ll/7“.  In  this  case  in 
the  transport  regions  will  be  comparable  to  or  larger  than  ambient.  We  shall 
see  below  however,  that  it  will  probably  not  be  possible  to  distinguish 
from  a  much  larger  component  R^  in  the  transport  regions. 
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Next  consider  fig.  D-2.  Results  for  R**  are  shown  on  fig.  D-2.  These  results  are 

analogous  to  those  for  R*  which  are  shown  on  fig.D-lln  that  they  are  horizontal 

profiles  and  R^J  at  z  -  90  km  =  ±  20  m.  Obviously  predosing  dominates 

in  the  region  PR.  Transport  dominates  the  regions  LTR  and  ATR.  Near  to  the 

region  IDR,  but  inside  the  regions  LTR  and  ATR,  we  expect,  on  the  basis  of 

the  same  arguments  that  are  applied  to  R^  in  the  transport  regions  that 
II  ^  -9  ,2 

R  =  4.5  x  10  watts/cm  sr  should  probably  be  close  to  the  actual  value 

II 

which  obtains  in  these  transport  regions.  The  R  in  the  transport  regions 
is  considerably  larger  than  R^  and  is  also  considerably  larger  than  the  maxi¬ 
mum  expected  ambient  horizontal  radiance  value. 


On  figures  D-3  &  D-4  we  show  vertical  cuts  at  x  =  2b,  0  and  -26  m  of  R^ 

and  R™.  These  show  that  we  can  expect  a  variation  of  the  order  0.5  to  1.5 
“9  2 

x  10  7  watts/cm  sr  on  the  interval  0>(z  -  90  km)  >  -  70  m.  It  may  be 
pose  ibis  to  detect  a  gradient  of  this  magnitude. 


Results  for  N0+(v)  Radiance:  We  expect  N0+(v)  will  be  strongly  quenched  by 

+  Mv)  4=  + 

N„.  +  NO  (v)  -»  N  +  NO  ( v- 1)  .  (3) 


From  Oordietz  et  al.  ("Theoretical  Study. . .Mechanisms, "  Mockba  ,  197b)  we  get 

"  3/2 

a  rate  constant  k^( l)  =  4.9  x  10  '  T  for  reactions  (2).  From  Tarr  and 

Archer  (DNA  5297F)  we  get 
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Figure  D-3  Estimates  R,_  for  the  Component  of  the  4.3pm  Radiance 
Vertical  Profiles  Which  are  due  to  the  Mechanism  I. 
These  vertical  cuts  are  for  x  =  -26,  0  and  26  m.  For 
z  <  90  km  there  is  essentially  no  x  variation  in 

(z).  The  large  variation  of  r|A°  in  the  transport 
region  where  z  -90  km  >  0  with  respect  to  R^D  in  the 
transport  region  where  z  -90  km  <  0  is  an  artifice  of 
the  approximate  method  which  is  used  in  the  calculation 
of  We  expect  10-^  watts/cnPsr  should  be  represen- 

t ive  of  the  actual  value  of  R*  which  should  obtain  in 


both  transport  regions. 


MECHANISM  II 


EXCITATION  OF  Nj 

N2  +  C°2  -  Nj  *  C°2 


PREDOSED  REGION  k  -  -H  m 


INSTANTANEOUS  DEPOSITION  a  -  • 


TRANSPORT  REGION  a  -  Mm 


TRANSPORT  REGION  a  -  -26,  0,  AND  4*6  * 

1+4/11/74 
~  BACKGROUND 


(«  -  90  km)  -  -50  m  -20  0  *20  40 

Figure  D-4  Estimates  for  the  Component  of  4.3ym  Vertical 

Radiance  Profiles  Which  are  Due  to  Mechanism  II, 
Excitation  of  N2;  N2  +  CO2  -*■  N2  +  CO2  (^3).  These 
vertical  cuts  are  for  x  -  -26,  0  and  26  m.  The  large 
variation  of  in  the  transport  region  associated 

with  z  -90  km  >  0  as  compared  to  R^J^  in  the  trans¬ 
port  region  where  z  -90  km  <  0  is  an  artifice  of 
the  method  by  which  the  approximate  solution  2 

is  calculated.  We  expect  rW  ^  4.5  x  10-^  watts/cm 
watts/cm^sr  should  be  near  the  actual  radiance  value 
in  the  transport  region. 


V 


kjOO/kjti) 


1.77  2 

1.27  3 

0.91  4 

0.59  5 

O.36  6 

0.22  7 

O.13  8 

0.07k  9 

0.043  10. 


We  also  get  from  DNA  3297F  an  efficiency  estimate  of  5  N0+  vibrational  quanta 
formed  per  ionization  event.  Equivalently,  £jjq+  *  *05  Is  the  fraction  of  the 
deposition  energy  converted  into  N0+  vibrational  excitation.  Finally,  from 
DNA  3297F  we  estimate  that  Just  the  N0+  bands  3  -*  2,  2  -»  1  and  1  -♦  0  may  be 
expected  to  contribute  radiation  which  is  overlapped  with  CO^  radiation  in 
the  4.1  to  4.4  um  region.  The  fraction  Fj-  of  the  total  N0+(v)  emission 
spectrum  contained  in  these  bands  is  F  —  ~  0.5  (DNA  3297F) . 


1  2  1 

From  Billingsleys  calculation  we  get  Aq  =  13-5/sec  and  that  A^  ■  2Aq,  we  assume 
3  1 

A^  =  3Aq  may  not  be  too  bad  an  approximation.  The  branching  ratio  for 

emission  rather  than  quenching  is 

3 

g  ~  _ 2 _ 

e/q  A^  +  k3(l)  [N2] 

and  with  k^(l)  ■  4.9  E-15  x  (220)^2  «  1.6  E-ll  cnP/sec  and  [N^]  »  10^'^^cm  ^ 
at  z  *  90  km  we  get 

Vq  ”  •°‘‘5 

So  finally  we  get 


111 


1  120  kv 

RN0+  “  >v  €N0+  F  3  DL  e/q 

R^q+  -  8  x  10-7  watts/cm2sr. 

This  Is  the  value  R^+  that  will  obtain  in  the  region  IDR.  Since  k^(l)[N^]  ~ 
WO/ nee  the  value  f  K^q+  will  rapidly  decay  in  the  region  such  that 
RN0+  <  tQr  x  <  '  l‘*-5  m.  Only  moleculELr  diffusion  will  contribute 

to  Rjjq+  the  transport,  regions,  we  have  not  as  yet  included  molecular 
diffusive  effects  into  our  calculation. 

Finally  we  emphasize  that  our  estimate  for  is  based  on  data  (e  +  and 

k^(v)/k^(l)  from  DNA  3?97F,  k^(l)  from  Gordietz  and  A^  ^  from  Billingsley) 
which  are  in  some  cases  unconfirmed.  Major  discrepancies  in  this  data  base 
could  induce  proportionate  discrepancy  in  our  estimate  for  R-^+.  Indeed, 
analysis  by  Herb  Mitchell  of  previous  rocket  borne  electron  gun  induced 
atmospheric  emissions  indicates  our  estimate  of  R^q+  which  is  based  on  the 
data  base  cited  above,  is  probably  too  large.  This  point  was  discussed  above 
in  somewhat  more  detail. 

In  our  study  we  did  not  consider  in  detail  quenching  due  to  N0+(v)  +  e  -♦  N  +  0 

but  if  we  take  a  =  3  x  10  cm  /sec  at  z  =  90  km  we  get  a[e]  =  5l*0/sec.  This 

is  comparable  to  k^(l)  [N„]  =  900/sec  so  apparently  recombinative  quenching  of 

N0+(v)  will  be  important  for  interpreting  the  EXCEDE  experimental  t.3  urn  data. 

*7  5 

Also,  we  expect  a(v)  may  actually  be  quite  a  bit  larger  than  3  x  10  cm  /sec, 
thus  increasing  the  importance  of  recombinative  quenching.  Recombinative 
quenching  is  not  important  for  the  auroral  case  however,  since  the  volume 
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ionizations  rate  is  down  by  about  six  orders  in  that  case,  hence  a[e]  is 
reduced  by  three  orders  of  magnitude. 


Summary  of  Results,  Discussion:  The  results  are  summarized  in  table  4-1. 

Clearly  R^q+  dominates  the  region  IDR  and  R^  dominates  the  transport  regions 

and  the  region  PR.  From  the  data  in  IDR  we  can  get  e  since  £„  +  *  R„,^+. 

NO  NO  NO 

From  data  in  the  region  PR  we  can  get  ,  the  efficiency  f  r  prompt  high 
energetic  election  coLliaion  induced  excitation  of  it  by  means  other  than  vv 
transfer  from  NO+(vJ  since  'c  R^.  Finally  we  may  get  s  >me  infor¬ 

mation  on  e  .  the  prompt  high  energetic  electron  collision  induced  efficiency 

for  prompt  eLectron  beam  excitation  of  CO  ,(v,)  from  data  in  the  transport 

I  II 

region  since  in  these  regions  R  /R  ir  0..1  '  based  on  the  off ic ienc ies  we  have 
assumed  In  trder  to  generate  these  preliminary  estimates,. 

Data  in  the  re*:  l  n  PR  will  help  t  •:  -nfirm  *hat  our  1  dimensi  snal  CO  urn 

time  dependent  radiati  'n  transport  -an  be  applied  with  confidence  4  the 
nuc  i  e'ir  cas.e  (rr-nrly  unif  rm  x-ray  patch}  .  Data  in  t  he  transp  ,rt  r>  gi  ms  can 
guide  us-  •'  that  we  tan  devel  p  a  numerically  effect  ive  '0  •*./;  am  5  dimen¬ 

sional  radiation  transp  -rt  facility  f >r  upplieatt  ci  t  4  the  nuclear  case  (highly 
structured  d  deposition  regi  ai) . 

Finally,  we  need  to  upgrade  ur  calculati  ns  it.  appli-atim  to  FXOFrK  as 
fa:  -ws : 

•  Generate  aa  iculnti  ns  f  >r  attitudes  sther  than  •»;  km. 

A 

•  Dev.'  1  p  capability  t.o  calculate  c,  rather  *  nan  the  less  ao  -usat <•  radial  ivi 

AD  ’ 

ilffu.  i  n  •  pps  >:imat.i  n  S._. 

AD  1 

and  pr  ve  t  be’  simil.ar  we  can  deve !  p  a  diftjsive  o.ppr  ximat  i  n 


•  If 


su  that  will  effectively  wipe  out  the  present  disparity  in  our  calcula¬ 
te  TT  T 

tions  for  R  ^  and  R*  in  the  region  LTR  as  compared  to  the  region  ATR. 

Such  an  approximation  S,^  may  be  necessary  to  adequately  account  for  struc¬ 
ture  in  the  t3  patch. 

Utilize  a  more  realistic  geometry  f>r  the  deposition. 

Account  for  molecular  diffusion. 


Finally  we  briefly  comment  on  what  we  expect  to  be  the  major  differences  in 
the  80  and  100  km  cases  In  comparison  with  the  90  km  case  that  we  have  treated 
here.  At  80  km  most  of  the  C0o  ^.3  u®  signal  we  see  will  be  dominated  by 
radiation  transport  effects,  the  region  IDR  will  be  just  7  m  in  length,  RN0+ 
will  be  significantly  depressed,  the  ambient  radiance  will  be  increasingly 
bothersome.  At  100  km  radiation  transport  effects  will  be  very  weak  as  will 
the  ambient,  radiance,  wil]  become  much  stronger. 


Appendix  E 

ESTIMATED  SYNTHETIC  SPECTRA 

Spectral  radiance  denoted  by  the  symbol  n  ®  1  to  6  on  the  figures  is  due 
to  lines  in  the  band  complex  n  as  defined  above.  The  symbol  n  ■  7  designates 
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Appendix  F:  Vufoils  Shown  at  the  1979  Spring 
Meeting  of  the  A.G.L'. 
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ANALYSIS  FOR  26  OCT  1978  ROCKETBORNE  MEASUREMENTS  OF 
4.3-pm  AURORA  AND  AIRGLOW  _ 


Table  I.  TICE  V.TAK  COg 

v.  HANDS 

y 

3 

00_(jv  )  -  (j) 

2  P 

Transition 

Isotope 

2 

0  0 

00  1  -  00  0 

636 

3 

0  0 

00  1  -  00  0 

628 

4 

0  0 

00  1  -  00  0 

627 

5 

Ol1!  -*  O^O 

626 

6 

Ol1!  -  0110 

636 

7 

Ol1!  -  Ol1© 

628 

8 

orli  -  oi^-o 

627 

9 

2  2 

02  1  -  02  0 

626 

10 

02°1  -  02°0 

626 

n 

10°1  -  10°0 

62  6 
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ROCKET  A  SCOUT 


JOOOO  4278A  BRIGHTNESS  (RAYLEIGHS) 

1  O'9  4.3  RADIANCE  (W/cm"2-*”1) 


.7 -pm  RADIOMETER 


MODELING  THE  10/26T78  UPLEG  DATA 


OH  COMPONENT 


TIME  BACK  FROM  PENETRATION  IN  MIN 


Appendix  G:  Calculation  of  Cell  Wheel  Interchange 
Weighting  Functions 


In  this  appendixwe  discuss  calculations  of  cell  interchange  weighting  func¬ 
tions  for  an  absorption  cell  wheel  filtered  radiometer  that  are 

analogous  tc  those  given  on  Fig.  5-9  in  the  article  "The  Infrared  Temperature 
Profile  Radiometer  (ITPR)  Experiment"  by  Smith  et  al. ,  which  appears  in  tne 
NASA  Nimbus  5  User's  Guide,  1972.  A  weighting  function  may  be  generated 
to  calculate  the  change  in  signal  R*^  that  is  observed  on  replacing  absorp¬ 
tion  cell  i  by  cell  j.  The  form  of  the  weighting  function  (o  N)  then  is 


A  ^ 

W  J  "  JV*"’  b  cp(&'  bx)SS  8Bmgm  e 

m  B 


(e  -  e 


(G-l' 


The  quantity  b  «  b  (z)  *  yT(z^) /T(z)  accounts  for  the  temperature  dependent 

altitude  variation  in  the  Doppler  line  width.  A  reference  altitude  z  ■  l«Otan 

is  utilized  in  the  calculation,  T(z^)  and  T(z)  are  the  atmospheric 

temperatures  at  z  and  z  respectively.  The  1/orentz  parameter  a  «  a(z)  ■ 
y  ■'  ■  fL 

(y0c  V  275  K/T(z)  p(z) ) /(vj)/b(z) )  accounts  for  pressure  and  temperature 

dependent  altitude  variation  in  the  collisional  contribution  to  the  line  shape. 

We  will  take  the  yQ  ■  .067  cm  ^atm  1  as  is  recommended  by  Ludvig  et.  al.,  NASA 

SP-5000,  1972.  The  quantity  p(z)  is  the  atmospheric  pressure  in  atm  at  z,  c 

is  the  speed  of  light  and  is  l/2  the  e  fold  Doppler  line  width  calculated 

for  the  temperature  T(z.).  The  quantity  cp(a(z) ,  b(z)x)  is  the  Voigt 

absorption  line  profile  at  z.  The  quantity  e^_  is  the  normalized  line 

Bm 

strength;  it  may  be  calculated  for  the  branches  B  ■  P,Q  and  R  by  formulas 

(9)  through  (12)  given  by  McClatchey  et.  al.,  on  page  1}  of  their  report 

v 

AFCRL-TR-73-0096.  The  sum  o  s_  »  1.  The  line  index  no.  m  is  defined  on 

d  Bm 

page  11  of  AFCRL-TR-73-0096 .  The  quantity  g^  «=  ®m(z)  the  fractional  pop¬ 
ulation  of  the  rotational  ground  state  corresponding  to  index  m,  i.e. 


g^z)  —  (2m+l)exp(-  he  Bm(m+l)/kT(z))/(vkT(z)/2hc  B  ~\fn' )  (C-2) 
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where  hc/k  —  1.439  cm/°k  and  the  C02  rotation  constant  B  —  0. 390218  cm"*. 

The  quantity  is  the  atmospheric  C02  optical  depth  above  z  at  x 

dimensionless  frequency  units  from  the  center  of  the  line  designated  by  Bm. 

The  quantity  x  -  (v^  -  v)/vD  where  is  the  Bm  line  center  frequency.  The 

quantity  w  -  l/2  for  the  626  001  -*  000  band.  The  quantity  t*,  is  given  by 

Bm  * 

00 

TBm(z>x)  "  0  J  dz'[C02]z*b(z‘’)Cp(a(z'),  b(z')x)  6Bmgm(2')  •  (G-3) 


The  quantities  t^x)  and  r^x)  are  the  optical  depths  across  the  ith  and 
Jth  cells  at  frequency  x  from  the  center  of  the  Bm  line.  They  are  given  by 


i  /  x  „i  .1,1  ,  i  N  i 

’Bm(x)  “  oN  b  cp(a  ,  b  xjs^ 


(G-4) 


where  N  is  the  C0o  column  density  across  the  ith  cell,  b  *  Vt(0/t 

i  ^  i  ^ 

where  T  is  the  temperature  of  the  ith  cell,  a  is  the  Lorentz  parameter  for 

the  ith  cell  which  is  determined  by  T*,  [CO,,]*  and  [B] 1  where  [(XL]*  and 
i  ^-2 

fB]  are  the  number  densities  of  C0o  and  a  possible  broadener  gas  B  in  the 
i  i  d 

ith  cell,  g  is  determined  by  T  via  equation  (G-2)  above. 


Ifce  change  in  nadir  radiance  R  which  occurs  when  cell  i  is  replaced 
by  cell  j  is  given  by 


%  RiJ  -  J*  hj  dN'  WiJ(cN#)Q(N')S(oN') 


(G-5) 


where  n(N')S(oN')  is  defined  in  the  paper  by  Kumer  and  James  (j.  Geophys. 
Rea.  J2,  638,  1974),  it  is  proportional  to  the  atmospheric  C02  026  001  -♦  000 
4.3  urn  volume  erission  rate.  Implementation  of  eqns.  (G-l)  through  (G-5) 
resulted  in  the  calculations  the  cell  Interchange  weighting  functions 

and  the  corresponding  cell  interchange  nadir  auroral  signals  that  are 

A 

shown  or  Figures  3-6  ard  3-7  respectively  in  the  main  text. 


Appendix  H:  Arc  Overfly  Nadir  Viewing  HIRIS  NESR  Calc. 


o  Av£  =  60  cm  3  cold  filter  width,  80%  transmission 
o  Resolution  Av=  0.5  cm  ^  achieved  by  stroke  of  Michelson 
o  Aperture  ^2.54  cm  (diameter) 

o  FOV  a  -  (2/(2347  cm-1/. 5  cm_1))3/?  -  2.07  x  10-2  i20.7«nrad 
m 

o  Fill  up  FOV  with  250°K  black  body  earth 

o  Bv(250°,  2347  cm-1)  =  3.742-12  (v3/ir) (eTO/250-l)"1  -  2.1  x  10-8  w/cm2 
str  cm~l  where  v  =  2347  cm-^  and  TO  ■  1.439v 

o  HIRIS  throughput  4  R  *  (w  x  (2.54/2)  x  (21  x  10_3/2) ] 2  x  0.4  - 
7  x  10  ^  (assume  0.4  optics  transmission) 

-20 

o  Photon  energy  hv  *  4.66  x  10  Joule 

o  Af  =  120  =  N (//  of  resolution  elements) 

o  Quantum  efficiency  n  =  0.2  at  4.3ym  for  A  :  S. 

s  I 

o  NER  *  (Bv  Av  $  Af/nhv)3^2(~r^)  «  4.5  x  10  ^  w/cm2  str  (PC  formula) 

H 

o  NESR  -  NER/Av  1  9  x  10~10  w/c.m2  str  cm-1 

This  NESR  is  achieved  in  1  sec.  of  integration  time.  Note  that  use  of  a  field 

widened  interferometer  would  decrease  the  NER  as  computed  above  by  the  ratio 

of  the  Improvement  in  affl.  Other  significant  improvements  by  use  of  the  field 

widened  interferometer  could  be  achieved  by  use  of  an  In:Sb  PV  detector 

rather  than  the  A  :S.  PC  detector, 
s  i 


184 


Appendix  I:  Contractually  Supported  Publication8 


Kumer,  J.  B.,  "Lateral  Structure  in  the  4.3pm  Arc",  EOS  58,  460,  1977 

Kumer,  J.  B.,  "Approximate  and  Exact  Technique  for  Multidimensional  Radiation 
Transport  in  a  Hane  Parallel  Atmosphere;  Application  to  the  4.3pm  Auroral 
Arc",  J.  Quant.  Spectrosc.  Radlat.  Transfer  19,  649,  1978 

Kumer,  J.  B. ,  A.  T.  Stair,  Jr.,  Ned  Whiller,  K.  D.  Baker  and  D.  J.  Baker, 
"Evidence  For  an  OH  W  Xv  (X^Cv^)  -*€02  +  hv  (4.3pm)  Mechanism  for  4.3pm 
Airglow",  J.  Geophys.  Res.  83,  4743,  1978 

Baker,  D.  J. ,  G.  D.  Frodsham,  J.  B.  Kumer,  A.  T.  Stair,  Jr.,  J.  C.  Ulvick, 
"Rocketbome  Measurements  of  Infrared  Aurora  and  Airglow  Emissions  at  4.3pm", 
EOS  60,  328,  1979 

Kumer,  J.  B.,  D.  J.  Baker,  G.  D.  F.  Frodsham,  A.  T.  Stair,  Jr.,  J.  C.  Ulvick, 
"Analysis  for  Measurements  of  4.3pm  Airglow  and  Aurora;  Evidence  for  a  Prompt 
4.3pm  Auroral  Component",  to  be  submitted  to  J.  Geophys.  Res.,  Oct.  1979 
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